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We present the universal theory of arbitrary, localized impurities in a confining paramagnetic state 
of two-dimensional antiferromagnets with global SU(2) spin symmetry. The energy gap of the host 
antiferromagnet to spin-1 excitations, A, is assumed to be significantly smaller than a typical nearest 
neighbor exchange. In the absence of impurities, it was argued in earlier work (Chubukov et al. Phys. 
Rev. B 49, 11919 (1994)) that the low-temperature quantum dynamics is universally and completely 
determined by the values of A and a spin-wave velocity c. Here we establish the remarkable fact 
that no additional parameters are necessary for an antiferromagnet with a dilute concentration of 
impurities, ni mp — each impurity is completely characterized by a integer/half-odd- integer valued 
spin, S, which measures the net uncompensated Berry phase due to spin precession in its vicinity. 
We compute the impurity-induced damping of the spin-1 collective mode of the antiferromagnet: 
the damping occurs on an energy scale F = ni mp (7ic) 2 /A, and we predict a universal, asymmetric 
lineshape for the collective mode peak. We argue that, under suitable conditions, our results apply 
unchanged (or in some cases, with minor modifications) to d-wave superconductors, and compare 
them to recent neutron scattering experiments on YBa2Cu3C>7 by Fong et al. (Phys. Rev. Lett. 
82, 1939 (1999)). We also describe the universal evolution of numerous measurable correlations as 
the host antiferromagnet undergoes a quantum phase transition to a Neel ordered state. 
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INTRODUCTION 



dinlcnsional doped antiferromag- 
in quantum many-body the- 
intciiesting quantum ground states 
itions appear possible, and the 
:perimental applications to the 
superconductors and other lay- 
mpounds. However the zoo of 
o the complexity of the prob- 
d quantitative theory has not 

a detailed quantum theory of 
of doped antiferromagnets. We 
which it is possible to neglect 
spin and charge degrees of free- 
of the spin excitations alone, 
in mind the following physical 
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systems, which are the focus of intense current experi- 
mental interest: _____ 

(A) Quasi-two dimensional 'spin gap' insulators-ffl-'- 
like SrCu203, CuGeC>3 or NaV^Os in which a small frac- 
tion of the magnetic transition metal ions (Cu or V) are 
replaced by non-magnetic ions by doping with Zn or Li. 
These insulators have a gap to both charge and spin ex- 
citations, but the spin gap is significantly smaller than 
the charge gap, validating a theory of the spin excitations 
alone. Moreover, additional low energy spin excitations 
are created in the vicinity of the dopant ions while charge 
fluctuations are strongly suppressed everywhere. 

(B) High temperature superconductors like YBa2Cu3C>7 
in which a small fractk^of-t^pLas. JacGa— replaced by 
non-magnetic Zn or L»__&_BBM_. A major 
fraction of the spectral weight of the spin excitiations 
near momentum Q = (tt, ^. resi d es , in a, spin- 1 'resonance 
peak' at energy A = 40meVBll_'E_l_3E_. We present here 
a theory for the changes in this peak due to the doping by 
non-magnetic ions. It may appear surprising that we are 
able to model this phenomenon by a theory of the spin 
excitations alone, but we shall describe how this is possi- 
ble: under suitable conditions, the coupling between the 
bosonic spin-1 mode responsible for the resonance peak, 
and the spin-1/2 fermionic, Bogoliubov quasiparticles of 
the d-wave superconductor is weak and, in a sense we will 
make precise, irrelevant. (The alternative case in which 
the coupling to the quasiparticles is relevant is discussed 
in Appendix |A| — here we find that qualitative features of 
the theory and the structure of all the scaling forms re- 
main unchanged, and there are only quantitative changes 
to scaling functions.) We will obtain a universal expres- 
sion for the energy scale over which the resonance peak is 
broadened by non-magnetic ions, with the result depend- 
ing only upon known bulk parameters. We also predict 
an unusual asymmetric lineshape for the broadened peak 
at low temperatures, and it would be interesting to test 
this in future experiments. 

Related questions have been addressed in some earlier 
work. The impurity-induced damping of bulk excitations 
of an antiferromagnet was addressed in Ref. however, 
their analysis was restricted to the doublet spin-wave ex- 
citations of an ordered Neel state well away from any 
quantum critical point; our focus here is on the triplet 
excitations in the paramagnetic phase. A weak-coupling 
analysis of the effect of Zn impurities on the spin fluctu- 
ation spectrum of a traditional BCS superconductor was 
considered in Ref. |2^. 

An elementary introduction and a summary of our re- 
sults, mainly for the uniform spin susceptibility, have ap- 
peared previouslyE_. 

Although our results apply to disparate physical sys- 
tems, they are unified by their reliance on a simple, new 
quantum field theory. We find it convenient to introduce 
this field theory at the outset, and will attempt to keep 
the discussion here self-contained and accessible to exper- 
imentalists. The discussion in this introduction will be 
divided into four subsections: Section I A will review the 



qua ntum field theory for the host antiferromagnet, Sec- 
tion I B will introduce the new quantum field theory of 



the impuri ty p roblem, our main results will be described 



in Section I C , while the outline of the remainder of the 



paper appears in Section [I] 



A. Host antiferromagnet 

We shall assume that the spin fluctuations in the host 
antiferromagnet, prior to doping by the Zn or Li ions, 
are described by the following 0(3)-symmetric quantum 
theory_3 of a three-component vector field (f> a (x, r) (a = 
1...3) 



Z' h = J 2?0 Q (x,T)exp(-S;) 



d 2 x I dr 



2 {(dr<p a ) 2 + cj(d x M 2 



^(fl_,*_) 2 + «tf_) + g(rf_) 



(1.1) 



Here repeated indices a are implicitly summed over, x = 
(xi,X2) is the two-dimensional spatial co-ordinate, r is 
imaginary time, 







n 



k B T> 



(1.2) 



and T is the temperature. The field cf> a (x,T) represents 
the instantaneous, local orientation and magnitude of the 
antiferromagnetic order parameter. So in applications to 
YBa2Cu3C>7, <p a represents the amplitude of spin fluc- 
tuations near wavevector Q = (7r,7r). More generally, 
(j) a can represent collinear spin fluctuations at any com- 
mensurate wavevector, apart from wavevector (0,0). So 
states with significant ferromagnetic fluctuations are ex- 
cluded, as are spiral states (like those on the triangular 
lattice antiferromagnet) which have non-collinear spin or- 
der. Incommensurate, but collinear, spin correlations can 
be treated by a simple extension of the theory we consider 
here - this is discussed briefly in Appendix pi 

The quantum dynamics of 4> a is realized by the second- 
order time derivative term in S^; we have rescaled the <j) a 
field to make the co-efficient of this term unity. Notice 
that there is no first-order Berry phase term, as is found 
in the path integral of a single spin: these are believed 
to efficiently average to zero on the scale of a few lattice 
spacings, because the orientation of the spins oscillates 
rapidly in the antiferromagnet. In one-dimensional anti- 
ferromagnets, this lattice scale cancellation is not quite 
perfect, and a topological '#-term' does survivecj; how- 
ever no such terms appear in two dimensions-the fate of 
the Berry phases in the hosL^—tiforromagnet has been 
discussed at length elsewhereooEj, and will not be im- 
portant for our purposes here. 

The spatial propagation of the spin excitations is as- 
sociated with the two spatial gradient terms, and we 
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have allowed for two distinct spin- wave velocities, Ci, 
c 2 , for propagation in the x\ and x-i directions. Such 
anisotropics are present in the coupled 'spin-ladder' sys- 
tems like SrCu 2 03. However, a simple redefinition of 
spatial co-ordinates 



ci 

x\ -> 1 — 

c 2 



1/2 



C2 

^2 -> 1 — 
Cl 



1/2 



£2, 



(1.3) 



allows us to scale away the anisotropy, and we shall as- 
sume this has been done in our subsequent discussion, 
unless explicitly stated otherwise. The resulting parti- 
tion function has the form 



5 b = 



where 



Z h = 

d d x / dr 
Jo 



J V<j) a (x, t) exp (-5 b ) 



~({d T cf> a ) 2 + C 2 (V x ct>a) 2 +S<f>l) 



"4! 



(cic 2 ) 



1/2 



(1.4) 



(1.5) 



We have also generalized the action to d spatial dimen- 
sions for future convenience. 

The most important property of 5b is that it|-exhibits 
a phase transition as a function of the coupling sc3. For s 
smaller than a critical value s c , the ground state has mag- 
netic Neel order and spin rotation invariance is broken 
because (4> a ) 0. For s > s c , spin rotation invariance is 
restored, and the ground state is a quantum paramagnet 
with a spin gap. The properties of this phase transition 
in 5b are well understood, and we review a few salient 
facts. The Neel state of 5 b is characterized by two spin 
stiffnesses, p s \, p s2 , which measure the energy cost to 
slow twists in the direction of the Neel order in the 1,2 
directions respectively (more precisely, a uniform twist 
by an angle 9 over a length L in the i direction costs 
energy (p S i/2)(9/L) 2 per unit volume; note that p S i has 
the dimension of energy in d = 2). These stiffnesses are 
related by 



(1.6) 



Psl _ c\ 

For the isotropic model 5b we have a single spin stiffness 



Ps = {PslPs2) 



1/2 



(1.7) 



As s approaches s c from below, the velocities ci and c 2 
remain constant, while 



Ps-iSc-s)^- 1 ^. 



(1.8) 



where v is a known exponent. The quantum paramagnet 
for s > s c is characterized by its spin gap A, and this 
vanishes as s approaches s c from above as 



As has been discussed at length elsewhereE30, pro- 
vided \s — s c \ and T are not too large, the energy scales 
A and p s , and the velocities ci, c 2 , are sufficient to com- 
pletely characterize the quantum dynamics of a d = 2 
antifcrromagnet; there is no need to have an additional 
parameter determining the strength of the quartic non- 
linearity go because it reaches a universal value, deter- 
mined by the zero of a renormalization group beta func- 
tion, near the critical point. The present paper will es- 
tablish the remarkable fact that no additional parameters 
are needed to characterize the spin dynamics in the pres- 
ence of a dilute concentration of impurities. As we shall 
specify more explicitly below, we only need to know the 
concentration of impurities, ni mp , and for each impurity a 
spin S, which is integer or half-odd-integer; the value of S 
can usually be determined, a priori, by simple arguments 
based upon gross features of the impurity configuration. 

Many simple, physically relevant, lattice antiferromag- 
nets can be explicitly shown to have a quantum phase 
transition described by 5b- One of the most trans- 
parent is the coupled-ladder antiferr.nni a^nf t .-ghich has 
been much studied in recent workEjO't2lEi}t3'LJ. We 
will consider such a model in Section 
lationship to 5b will emerge naturally, 



III A 



and 



the re- 
ive will 



A~ (s 



(1.9) 



also be able to relate parameters in 5b to micro- 
scopic couplings—Other models include the double-layer 
antifcrromagnetaO't3EjEi-|the antifcrromagnet on the 
1/5-depleted square lattices, and .the square lattice an- 
tifcrromagnet with frustratiorOcd'OCj: our results will 
also apply to such antiferromagnets. 

The reader might wonder why we are placing an em- 
phasis on the quantum phase transition in 5b, while most 
experiments on Zn or Li doping have taken place in an an- 
tifcrromagnet or superconductor with a well-defined spin 
gap. The reason is that the spin gap, A, is often sig- 
nificantly smaller than the microscopic exchange interac- 
tions, J. Under such circumstances, it is a poor approx- 
imation to model the spin gap phase in terms on tightly 
bound singlet pairs of electrons. Rather, the smallness of 
the gap indicates that there is an appreciable spin corre- 
lation length and significant resonance between different 
pairings of singlet bonds. The approach we advocate to 
describe such a fluctuating singlet state is to find a quan- 
tum critical point to an ordered Neel state somewhere in 
parameter space, and to expand away from it back into 
the spin gap phase. It is not necessary that the quan- 
tum critical point be experimentally accessible for such 
an approach to be valid: all that is required is that it be 
near a physically accessible point. The value of expand- 
ing away from the quantum critical point is that detailed 
dynamical information can be obtained in a controlled 
manner, dependent only upon a few known parameters. 



B. Quantum impurities 

We now introduce arbitrary localized impurities at a 
set of random locations {r} with a density 7ii mp per unit 
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volume. We place essentially no restrictions on the man- 
ner in which each impurity deforms the host antiferro- 
magnet, provided the deformations are localized, immo- 
bile, well separated from each other, and preserve global 
SU(2) spin symmetry. As we argue below, each impurity 
will be characterized only by a single integcr/half-odd- 
integer, S r , which measures the net imbalance of Berry 
phase in the vicinity of the impurity; all other character- 
istics of the impurity will be shown to be strongly irrel- 
evant. We mentioned the essentially complete cancella- 
tions of Berry phases in the host antiferromagnet earlier 
in Section I A; this cancellation can be disrupted in the 
presence of an impurityC-JEJ, and it is only this disrup- 
tion which will be important in the low energy limit. For 
instance, a single Zn or Li impurity introduces one addi- 
tional spin on one sublattice than on the other, and so in a 
host spin- 1/2 antiferromagnet such an impurity will have 
S r = 1/2. As another example, consider a ferromagnetic 
rung bond in a spin-1/2 coupled ladder antiferromagnet: 
the two spins on the ends of the rung will be oriented 
parallel to each other and so their Berry phases will add 
rather than cancel, leading to S r = 1. 

The above criteria determining S r are qualitative, but 
we can make them more precise. Place the host antiferro- 
magnet in the quantum paramagnetic phase (s > s c ) and 
measure the response to a uniform magnetic field. In the 
absence of impurities, the total linear susceptibility of the 
antiferromagnet, x-, can be written as \ — (3/is) 2 'Axh, 
where A is the total area of the antiferromagnet, g is the 
gyromagnetic ratio, [Ib is the Bohr magneton, and Xb is 
the response per unit area. As T — > 0, the response is ex- 
ponentially small due to the presence of a spin gap, A; a 
simple argument summing over a dilute gas of thermally 
excited spin-1 particles shows thatEj 



Xb(T - 0) 



A 



nn c z 



-A/k B T 



(1.10) 



Now let us add a very small concentration of impurities. 
The response x is modified to 



X= {gVB) 2 [Axb + Ximp\ 



(1.11) 



Provided the impurities are dilute enough that they do 
not interact with each other, then Ximp must take the 
following form at low T: 



1) 



3k B T 



(1.12) 



We will take the values of S r which appear in such an 
observation of Ximp as the definitions of impurity spins 
S r at the site r: they will necessarily be integers/half- 
odd-integers when the host antiferromagnet has a spin 
gap. It is important to note that for any fixed n- lmp , 
there is an exponentially small temperature below which 
the impurities do interact with each other, and ( 1.12| ) is 
no longer valid: we have assumed that T is above such 
a temperature in determining the S r . However, once the 



S r have been obtained, our quantum theory below will 
also apply in the interacting impurity regime. 

We now describe the modifications to the bulk action 
5b by the quantum impurities with spi n S r at sites {r}. 
The divergent Curie susceptibility in ( 1.12p is a conse- 
quence of the free rotation of the impurity spin in spin 
space. Let us denote the instantaneous orientation of 
the impurity spin at r by the unit vector n ra (r), where 
a = 1 ... 3 and J2 a n ra( T ) = 1- Then the field theory of 
the quantum impurity problem is 

Z = J Vcj) a (x, T)Vn ra {r)8 {n 2 ra - l) exp (-<S b - Simp) 
Smp =z2 dr 



Of /Iq. 1 Thy J " 



dr 



lorSAa (x = r, r)n ra (r) 



(1.13) 



The first term in Si mp is uncompensated Berry phase 
of the impurit y at site r: the spin S r is precisely that 
appearing, in (1.12) and A a (n) is a 'Dirac monopole' 
functioned which satisfies 



(1.14) 



The coupling between the impurity spin and the host 
spin orientation is 7or and this can be either positive 
or negative-this depends upon the sublattice location of 
the impurity. The main purpose of the remainder of this 
paper is to describe the properties of Z in its different 
physical regimes. 



Magnetic 
Field H 



I t \ t I 
t I t I t 



♦ tot* 
t * 1 1 1 

1 1 i 1 1 



FIG. 1: From Julien et al. 



(Ref. 0). 



Polarization of the 



Cu spins around the central Zn impurity (the-icircle) in the 
presence of an external field H, as measuredlij by NMR on 
YBa2Cu3Cv7- The size of each arrow represents the value of 
the moment. Note that there is no spin directly on the Zn 
site, but it is expected that a net moment of S = 1/2 centered 
at it at T = 0. 

It is worth emphasizing that <Si mp is to be viewed as 
a long-wavelength theory for the impurity at site r, and 
requires some interpretation when applied to lattice mod- 
els. In particular when a non-magnetic Zn or Li impurity 
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replaces a spin-1/2 Cu ion, there is clearly no spin degree 
of freedom directly at the impurity site. Instead the spin 
polarization is distributed in a staggered arrangement on 
the neighboring Cu ions, so that the total moment is ex- 
pected to be S r = l/i? at f = 0; this has been studied in 
elegant experimentsoEia, and we reproduce the results 
of one of these observations in Fig. |l]. The unit vector 
n ra then measures the instantaneous orientation of the 
collective spin polarization which is cent ered on a non- 
magnetic Zn site, and the Berry phase in (1.13) is the net 
uncompensated contribution obtained by summing over 
the Berry phases of all the spins on the neighboring Cu 
ions. The spatial correlations of the spins on the Cu ions 
are encoded, in our theory, in those of the 4> a field. 

Som e crucial properties of the partition function Z in 
(1.13) follow from simple power-counting arguments. We 
will perform these here at tree level, and higher loop cor- 
rections do not lead to corrections which modify the main 
conclusions. We focus here on the properties of a single 
impurity: the consequences of interactions between the 
impurities will be discussed later. We examine the behav- 
ior of <Sb+iSi m p under a rescaling transformation x — > x/b 
and t — ► t Jb centered at the impurity (the dynamic crit- 
ical exponent, z, is unity). The scaling dimensions of <fi a 
and n a are fixed by demanding invariance of the time 
derivative terms: this leads to 

dim[</> Q ] =(d- l)/2 ; dim[n Q ] = 0. (1.15) 

These can be used to determine the dimension of 7o r : 

dim[ 70r ] = (3 - d)/2. (1.16) 

So 7or is a relevant perturbation about the decoupled 
fixed point in d jr 2 . F ollowing results of earlier work on 
related modelsc§e3E3, we find that in an expansion in 



e = 3 — d, 



(1.17) 



7or reaches one of two fixed point values ±7* under the 
renormalization group (RG) flow (the value of 7* depends 
upon S r , but we will not denote this explicitly). So the 
low energy properties of Z are actually independent of 
the particular bare couplings 7o r , and are controlled in- 
stead by 7*. The memory of the initial values of 70 is 
retained in the sign of the fixed point value: at the fixed 
point we have 7o r = 0>7* where a r = ±1, and the sign 
of 0> is the same as that of the initial sign of 70, i.e., 
7or does not change sign under the RG flow. The initial 
sign of 7o r is determined by sublattice location of the 
impurity in the host lattice antiferromagnet, and so both 
signs are equally probable. 

Similar arguments can be used to show that there are 
no other relevant couplings between the impurity and 
the host antiferromagnet, consistent with the required 
symmetries. The simplest possible new coupling is 



which represents variation in the strength of the exchange 
constants in the vicinity of the impurity. Power counting 
shows that 



dim[Cr] = 2-d. 



(1.19) 



So the ( r are strongly irrelevant for small e (this is all 
that is needed, as the tree level arguments are valid only 
for small e). We can also consider the coupling of the 
impurity spin to the local uniform (ferromagnetic) mag- 
netization, density L a . For the action 5b, the latter is 
given byE3 

L a = iea/3 7 0/3<9 r </> 7 ; (1-20) 
and the coupling to the impurity spin takes the form 



J2jdrCrSrL a{x = r,r)nr a (r). 



(1.21) 



E 



dr( r (j) 2 a {x = r, r), 



(1.18) 



Because L a is the density of a conserved charge 

dim[i Q ] = d, (1.22) 

and so 

dim[Cr] = 1 - d, (1.23) 

and ( r is also irrelevant. Similar arguments apply to all 
other possible couplings between the host and the impu- 
rities consistent with global symmetries. 

(It is worth mentioning, parenthetically, that if we re- 
lax the constraint of global SU(2) symmetry, and allow 
for spin anisotropy, then relevant couplings are possible 
at the impurity site. The most important of these are 
a local field ^ J drw r n rz {T), or a single-ion anisotropy 
J2 r I dTw r n 2 z (r): the w r , w r are easily seen to be rele- 
vant with scaling dimension 1 . Now there is one relevant 
coupling in the bulk (s), and one or more on the impurity 
(w r , Wr), and even the single-impurity problem has the 
topology of a multicritical phase diagram: separate phase 
transitions are possible at the impurity, the bulk, or both- 
and the full technology of boundary critical phenomenaEJ 
has to be used. For S r = 1/2 the single- ion anisotropy is 
inoperative, and if there is also no local field, we have to 
consider exchange anisotropics in the 7o r , and these can 
also be relevantc3.) 

We have now assembled all the ingredients necessary 
to establish our earlier assertion that no additional pa- 
rameters are needed to describe the dynamics of the an- 
tiferromagnet in the presence of impurities. For the case 
of a single impurity (or a non-extensive density of impu- 
rities) our arguments show that there is only one relevant 
parameters at the quantum critical point — the 'mass' s; 
its scale is set by the spin gap, A, for s > s c and by the 
spin stiffness, p s , for s < s c . All other couplings, whether 
they are bulk (like go) or associated with the impurities 
(like 7or) either flow to a fixed point value or are strongly 
irrelevant. We do need a parameter to set the relative 
scales of space and time, and this is the velocity c. Our 
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arguments also allow us to obtain important information 
for the case where the density of impurities, n; mp , is fi- 
nite. Now n; mp itself is a relevant perturbation of the 
critical point; however, we can use the fact that there is 
no new energy scale which characterizes the coupling of a 
single impurity to the bulk antiferromagnet to conclude 
that the value of n; mp alone is a complete measure of the 
strength of this perturbation. 

Let us restate the above important conclusion a lit- 
tle more explicitly. Assume we are doping the spin gap 
paramagnet with s > s c and a spin gap, A (very similar 
arguments apply also for s < s c with p s replacing A). For 
the case of a single impurity, the dynamics of the host 
antiferromagnet in the vicinity of the impurity is deter- 
mined entirely by A and c with no free parameters. A 
finite density of impurities, ni mp , is a relevant perturba- 
tion, but its impact is determined entirely by the single 
energy scale, F, that can be constructed out of the above 
parameters, and that is also linear in n; mp : 



r 



a P (hc)° 



(1.24) 



In particular, all dynamic and thermodynamic proper- 
ties of the host antiferromagnet are modified only by a 
universal dependence on the ratio T/A; all corrections 
to these (including those dependent on the magnitude of 
the bare coupling between the impurity and host anti- 
ferromagnet) are suppressed by factors of A/ J, and only 
the latter are non-universal. The leading universal mod- 
ification depends also on the statistics of the values of 
S r and oy: in most physical applications we simply have 
S r = S independent of r, while the o~ r are independently 
distributed and take the values ±1 with equal probabil- 
ity. 

The above universal dependence on T/A does not pre- 
clude the possibility that there could be a phase transi- 
tion in the ground state as a function of T/A: a magneti- 
cally ordered state could appear for infinitesimal doping, 
i.e., at T/A = + , or at some finite critical value of 
T/A — these two possible phase diagrams are shown in 
Fig ||. Determination of the phase diagram requires solu- 
tion of the problem of interacting quantum impurities; in 
the paramagnetic phase (s > s c ), an effective action for 
interactions can be obtained by integrating out the <p a 
fields from Sb + Simp • Apart from the same Berry phase 
terms in iSi mp , the resulting action contains the term 



- Jdrjdr' S r S r io- r o- r i D(r- 



r')n ra (T)n r , a (T'), 



(1.25) 

where the universal interaction D is proportional to the 
4> a propagator and depends only on the fixed point value 
7* (which is an implicit function of S r ) — D decays ex- 
ponentially for large |r — r'\ or |r — r'|; there are ad- 
ditional multi-spin interactions between the impurities, 
wh ich h ave not been displayed. A notable propertyErCJ 
of (1.25) is that the signs of thq interactions have the dis- 
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Quantum 
paramagnet 



(b) 



FIG. 2: Two possible phase diagrams of the quantum field 
theory Sb + Sim P for a dilute concentration of impurities ni mp . 
For s > s c , all low energy properties depend only upon the 
dimensionless ratio F/A, where A is the spin g ap of the un- 
doped host antiferromagnet, and F is defined in (1.24). In (a) 
there is a phase transition to an ordered state at T/A = + , 



and a quantum paramagnet is present only for ni n 



and 



S > S c 



order of the "Mattis model 



"51 



i.e. 



there is no frustra- 



in (b) the transition to the ordered st ate occu rs at 
some finite universal value of F/A, and so by (1.9,1. 24 ) the 
phase boundary obeys 7ii mp ~ (s — s c )" d . For s < s c the prop- 
erties are universal functions of the analogous dimensionless 
ratio ni m p (Tlc / ps) 2 (where p s is the stiffness of the host anti- 
ferromagnet); we have not shown a third possible phase dia- 
gram in which the phase boundary is present for s < s c with 
limp ~ (s c — s)" d . Right at s = s c , the spacetime correla- 
tors are universal function of the dimensionless combinations 
xn im!p an d crn imp with no arbitrary scale factors. For any lat- 
tice antiferromagnet, the long-range order will disappear at 
sufficiently large ni mp when the concentration of impurities 
is near the percolation threshold — this phase boundary has 
not been shown in the figure, as it is not a property of the 
continuum theory <Sb + <Si mp . This-Jatter transition was con- 
sidered by Chen and Castro NetoEj, along with the behavior 
as a function of nj mp for s < s c ; however, they did not include 
the Berry phases,r-and these are expected to be important in 
the critical regionCJ. 



tion, and the product of interactions around any closed 
loop has a ferromagnetic sign. However, unlike the clas- 
sical spin glass, the fluctuating signs of the exchange in- 
teraction cannot be gauged away in this quantum Mattis 
model, as the transformation n a — > —n a is not a sym- 
metry of the Berry phase term (equivalently, it does not 
preserve the spin commutation relations). Nevertheless, 
the absence of frustration in such a modeLias been used 
by Nagaosa et alE2 and Imada and Iinoca to argue in 
favor of the phase diagram in Fig |2]a. 



7 



We reiterate that, irrespective of the correct phase di- 
agram, all low energy properties depend only the value 
of the dimensionless ratio T/A, which is a measure of 
the strength of the relevant perturbation of a finite con- 
centration of impurities. It is remarkable that such a 
measure is independent of the magnitude of the exchange 
coupling between the impurities and the antifcrromagnet. 
We also note that-some of the earlier scaling arguments 
of Imada and IincO are contained in such an assertion — 
it implies their identification of the crossover exponent of 
the impurity density as vd. 



1. Application to d-wave superconductors 

While it is reasonably transparent that the action 
Sh + <5i mp applies to the Zn doping of insulating spin gap 
compounds like SrCu203, the applicability to Zn doping 
of a good d-wave superconductor like YBa 2 Cu307 has 
not yet been established. Actually the similarity of Zn 
doping in a cuprate superconductor to that of spin gap 
compounds was already noted by the experimentalists in 
Ref 13— here we will sharpen their proposal. 

There is a great deal of convincing experimental ev- 
idence that each Zn impurity has an antiferromagnetic 
polarization of Cu ions in its vicinity, and the net mo- 
ment of this polarization is expected to be S r = 1/2 at 
T = 0; at moderate temperature this moment precesses 
freely giving rise to a Curie susceptibility, while there 
is evidence for spin freezing at very low temperatures, 
presumably from inter-impurity interactions (see Fig. [l], 
Refs. HJ|,|l^,|lJ and references therein). The Berry phase 
term in 5; mp then encodes the dynamics of the collective 
precession of the spins on the Cu ions near the Zn impu- 
rity at site r (recall that there is no spin on the Zn ion 
itself). 

The impurity spin will couple to spin excitations in 
the host superconductor. The most important of these is 
the spin-1 cnUectiyf .mode which gives rise to the 'reso- 
nance peak'OlHllj'EiEy of the cuprate superconductors. 
In our approaches this collective mode is represented by 
the oscillation of the field <j> a about (j) a = controlled 
by the action Sh- So our picture (Afche resonance peak 
is similar in spirit to computationsEj'LJLj which identify 
it with a S = 1 bound state in the particle-hole chan- 
nel in RPA-like theories of the dynamic spin response of 
a d-wave superconductor. However, such theories use a 
weak-coupling BCS picture of the electron spin correla- 
tions and also neglect the non-linear self-interactions of 
the collective mode. In contrast, in our approach the un- 
derlying spin correlations are characterized by 5b , which 
assumes that the superconductor is close to an insulating 
state with magnetic long-range order (and possibly also 
charge order): this has the advantage of allowing a sys- 
tematic treatment of the strongly relevant quartic non- 
linearity in the <p a field, and these non-linear effects are 
crucial to the universal nature oLour results. We also 
mention the approach of ZhangLJ which assumes that 



3-component 4> a is part of 5-component "superspin" — 
unlike Zhang, our theory does not appeal to any higher 
(approximate) symmetry group in the superspin action, 
beyond that expected from spin rotation invariance. 

Spin is also carried by the fermionic, spin-1/2 Bo- 
goliubov quasiparticles in a d-wave superconductor, and 
these will couple to the impurity spin. These quasiparti- 
cles have vanishing energy at four points in the Brillouin 
zone - (±-KT, ±K). As is well known, a gradient expansion 
the of low energy fermionic excitations in the vicinity of 
these points yields an action that can be expressed in 
terms of 4 species of anisotropic Dirac fermions. We do 
not wish to enter into the specific details of this action 
here, as only some gross features will be adequate for our 
basic argument. Let us represent these fermions schemat- 
ically by the Nambu spinors then the action has the 
form 

S,j, = J d 2 xdr [tft0 T tf + tf+Va*] • (1.26) 

We have omitted all coupling constants and matrices in 
the Nambu and spin spaces, as we do not need to know 
their structure here. The important point is that this 
action implies the dimension 

dim[#] = 1 (1.27) 

under scaling transformations. 

Now let us couple Sfr to the degrees of freedom in Sh + 

^imp • 

There will be a bulk cubic coupling of the form ^ a $$ 
or ifia^'i' only if it is permitted by the momentum con- 
servation in the host antiferromagnet. As the 4> a repre- 
sent spin fluctuations at the wavevector Q, such a cubic 
coupling is permitted only if Q equals the sum or differ- 
ence of two of (±K, ±K) {i.e. if Q = (2K, 2K) or not). 
For simplicity, in the body of the paper, we will assume 
in this is not the case, i.e. we assume Q ^ (2K,2K). 
Indeed, in YBa2Cu307, Q = (ir, n) and it appears that 
K 7^ 7r/2, and so a cubic term is forbidden. However, K 
is not too far from x/2, and so the effect of a cubic term 
may be manifest at higher energies. In Appendix |A| we 
consider the limiting case Q = (2K,2K), which permits 
a cubic term in the host d-wave superconductor at the 
lowest energies. A renormalization group analysis of such 
a host theory was presented recently by Balents et al^B, 
and Appendix ^ combines their results with the quan- 
tum impurity theory of the present paper: we find that 
the scaling structure is essentially identical to that in the 
body of the paper, and there are only simple quantitative 
modifications to the fixed-point values of the couplings. 

For now, we assume that momentum conservation pro- 
hibits a bulk coupling between the bosonic (4> a ) and 
fermionic (\I/) carriers of spinEll. However these degrees of 
freedom can still interact via their separate couplings to 
the impurity spins. The coupling of the superconducting 
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quasiparticles to the impurities at sites {r} is of the form 

S*im P = ^ / dT[J r 5 r n ra (T)^(x = r,T)*(x = r,r) 

+J' r S r n ra {T)^{x = r, r)*(x = r,r) + H.c. 
+K.* t (a; = r, r)*(x = r, r) 

= r,r)*(x = r,r) +H.c.]. (1.28) 

Again, we have omitted matrices in spin and Nambu 
space; the J r , J' r are exchange couplings between the 
quasiparticles and the Cu spins in the vicinity of the 
Zn impurity, and the V r , VI represent potential scatter- 
ing terms from the non-magnetic Zn ion itself. Mod- 
els closely related to + S^ imp have been the sub- 

jcct ^l^^l^rted f studics in thc past few 

yearSSSBElOOSBEl and a number of interesting 
results have been obtained. However, all of these earlier 
works have not included a coupling between the impurity 
spin and a collective, spin-1, bosonic mode like <p a . One 
of the central assertions of this paper is that (provided 
the spin gap, A, is not too large) such a coupling (as 
in iSimp) is of paramount importance for thc low energy 
spin dynamics, while the coupling to the superconducting 
quasiparticles (as in iS*i mp ) has weaker effects. 

(However, when one is explicitly interested in quasi- 
particle properties, as in tunneling experiments, it is cer- 
tainly necessary to include 5*.+ iS*i mp ; in STM experi- 
ments with atomic resolutionEEl, the quasiparticle tunnel- 
ing can be observed directly at the Zn site, and there the 
potential scattering terms V r ,V^ will be especially impor- 
tant as there is no magnetic moment (see Fig [l]). As we 
noted below (1.14), we are considering a long- wavelength 
theory of the spin dynamics, and careful interpretation 
is required for lattice scale effects.) 

T he ar gument behind our assertion is quite simple. Us- 
ing (1.27) we see that 



dim[J r , J' r ,VrX\ = -1. 



(1.29) 



Therefore, while the couplings, 7o r , between the impu- 
rities and the 4> a had a positive scaling dimension in 



(1.16), those between the fermionic quasiparticles and the 
impurities have a negative dimension and are strongly 
irrelevant. We are therefore justified in neglecting the 
fermionic quasiparticles in our discussion of the impurity 
spin dynamics. 

onclusion is also supported by 
of models related to <S* + 5*im P 
in the context of "Kondo problems with a pseudogap in 
the fermionic density of the states": for the case where 
the fermionic density of states vanishes linearly at the 
Fermi level (as is the case for Sy), there is no Kondo 
screening for small and moderate J r , J' r values, and the 
impurity spin is essentially static. Of course, the present 
renormalization group argument cannot rule out the 
possibility of new physics, associated with the fermionic 
Kondo effect, appearing at very large J r , J' r : we shall 
assume that the bare couplings are in a regime such that 



t his h as not happened, and the scaling dimensions in 
( 1.29 ) continue to apply. 



C. Results 

We will state our mai n results for the case of a single 
impurity in Section I C 1 , and those for a finite density in 
Section II C 2. 



1. Single Impurity 

We first consider a single impurity at the origin of co- 
ordinates x = with 



S r =a = S 



(1.30) 



An important measurable resp onse function is the im- 
purity susceptibility defined in (1.11). Its properties fol- 



low simply from a naive application of the scaling ideas 
we have presented above, and are summarized in Fig [jj. 
Knowledgeable readers may be surprised by such an as- 
sertion. In the intensively .studied multichannel Kondo 
quantum impurity problemo naive scaling actually fails: 
even though the low energy physics is controlled by a 
finite,-epupling quantum critical point, a 'compensation' 
effectL3 causes most thermodynamic response functions 
to vanish in the naive scaling limit, and the leading low 
temperature behavioE-4s-exposed only upon considering 
corrections to scaling 7 ttI. Fortunately, our problem is 
simpler — the bosonic 4> a excitations which 'screen' the 
impurity are themselves controlled by a non-trivial inter- 
acting quantum field theory 5b, and it is not possible to 
'gauge away' the effect of an external magnetic field on 
them. 

Armed with this reassuring knowledge, we need only 
know that Ximp has the dimensions of inverse energy to 
deduce its critical properties. For s > s c we of course 



have the Curie response in (1.12), which is, not coinci- 
dentally, also consistent with the scaling requirements. 
At the critical coupling s = s c , we continue to have a 
Curie-like response (because now fc^T is the only avail- 
able energy scale) but can no-longer require that the ef- 
fective moment be quantized: 



Xu 



Ci 

k B T 



T > 



(1.31) 



The universal number C\ is almost certainly irrational, 
and we will estimate its value in Sections II B and [II D ; 



our m ost reliable estimate for 5 = 1/2 and d = 2 is in 
( |3.36 ), with Cf TOC = 1/4. Turning finally to s < s c , the 
T = response is now anisotropic because of the Neel 
order in the ground state. We consider here only the 
response transverse to the N eel ord er, x±i at T = 0: 
further details are in Section II C 2 . The same scaling 
arguments now imply 



Ximp± 



c. 



; T = 0, d = 2, s<s c 



(1.32) 
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Xh = C 2 k B T/(h 2 c 2 ) / 



Ximp = c l/( k BT) | / 



Xb± = Ps/(Ji 2 <?) 
XimpT = C 3/Ps 



\ 

\ / 
\ / 
\ / 



&-A/k B T 
ttH c 2 

Ximp = 
S(S + 1)/(3* S T) 



FIG. 3: Summary of the results for the bulk and impurity sus- 
ceptibilities of 5b + <Simp in d = 2 with a single impurity at the 
origin of co-ordinates with S r = S. The dashed lines represent 
crossover boundaries. The crossover for s > s c occurs when 
A ~ ksT, and is between the high-temperature 'quantum 
critical' regime and the low-temperature quantum paramag- 
net described by a dilute gas of spin-1 (f> a quanta. Similarly, 
the crossover for s < s c occurs when p s ~ fcsT, and the low- 
temperature regime is a locally ordered Neel state whose long 
range order is disrupted by long- wavelength, classical, spin- 
wave fluctuations. The constants Ci_3 are universal numbers, 
insensitive to microscopic details. The constant C2 was intro- 
duced in Ref and is determined solely by the bulk theory 
Sh', quantum Monte Carlo results were obtained in Ref |3!j. 
The constants Ci and C3 are determined by Sh + <Si mp , and 
depend only on the integer/half-odd- integer valued S. The 
constant Ci defines the effective spin at the quantum-critical 
point by C\ = S e «(S e g + l)/3, and S e g is neither an inte- 
ger nor a half-odd-integer. On the ordered side, s < s c , in 
the absence of even an infinitesimal spin anisotropy, the mea- 
sured-spin susceptibility will take the rotationally averaged 
valueEa Ximp = S 2 /(3fcsT) + (2/3)xim P ±; we have no rea- 
son to expect non-monotonic behavior of the susceptibility 
as a function of s, and so this result suggests the bounds 
S 2 /3 < Ci < S(S + l)/3 in d = 2. 



Again C3 is an irrational universal number whose value 
will be estimated later. 

Having described the response to a uniform, global 
magnetic field, let us consider the responses to local 
probes in the vicinity of the impurity. Such results will 
apply to NMR and tunneling experiments. In consider- 
ing such response function it is useful to use the language 
of boundary conformal field theorjO: we are considering 
here a bulk (d+ l)-dimensional conformal field theory, Sh 
at s = s c , which contains a one-dimensional 'boundary' 
degree of freedom at x = 0. Correlation functions near 
the boundary will be controlled by boundary scaling di- 
mensions which are distinct from the scaling dimensions 



of the bulk theory we have considered so far. The most 
important of these controls the long-time decay of the 
impurity spin: 



(n a (T)n a (0)) 



1 



T = 0, t -> 00. (1.33) 



(Here, and in the remainder of the paper we are indulging 
in a slight abuse of terminology by identifying these as 
the correlations of the "impurity spin"; as in Fig [j] it 
is often the case that there is no spin on the impurity 
site — in such cases we are referring to spin correlations 
at sites very close to the impurity.) The quantity rf is 
an anomalous boundary exponent: we will compute its 
value in the expansion in e in Section II A. Clearly, (1.33) 
implies the scaling dimension 



dim[n Q ] = r//2, 



(1.34) 



which corrects the tree-level result in ( 1.15| ). (This is a 
good point to mention that loop corrections also m odify 
the scaling dimension of the bulk field <p a 



1.15) to 



6J = (d - 1 + 77)72, 



(1.35) 



where rj (like, v) is a known exponent which is a property 
of 5b aloneEj.) 

For s 7^ s c , there is a finite remnant moment in the 
boundary spin correlations. For s < s c , this is simply 
a consequence of the bulk Neel order also breaking the 
spin rotation symmetry on the boundary too. However, 
for s > s c , this is a somewhat more non-trivial quan- 
tum effect: the presence of the spin gap in the bulk, 
means that the boundary excitations, associated with the 
Berry phase term in 5i mp are confined to the impurity 
and maint ain a permanent static moment. So we may 
generalize (1.33) to 



lim (n a (r)n a (0)) = mf 



2 

imp 



T = 0. 



(1.36) 



The impurity moment, mi mp , (the remarks below ( 1.33| ) 
on abuse of terminology apply to the "impurity moment" 
too) behaves like 



m,. 



(1.37) 



a consequence of the scaling dimension (1.34). 

We have computed a large number of correlation func- 
tions describing the spatial and temporal evolution of the 
spin correlations in the vicinity of the impurity. We will 
leave a detailed discussion of these to the body of the 
paper, but note here some simple arguments which allow 
deduction of important qualitative features. The exper- 
imental probes mentioned earlier can follow the spatial 
evolution of the correlators of the bulk antifcrromagnetic 
order parameter field (f> a , and also that of the uniform 
(ferromagnetic) magnetization density L a . As the spa- 
tial arguments of 4>a(x, r) and L a (x, r) approach the im- 
purity, their critical correlations must mutate onto those 
of the boundary degrees of freedom. This transformation 
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is encoded in the statements of the operator product ex- 
pansion 

v a i \ n «( T ) 
lrm^r) ~ lxl(d _ 1+v _ rn/2 



hmi Q (^r) ~ j-^. 



(1.38) 



where the p owers of x follow immediately from the s caling 
dimensions (|l . 34| , p~3^ , []~2^) . The results fll.33|) and ( |l.3g| ) 



can be combined with scaling arguments to determine the 
important qualitative features of the spatial, temporal, 
and temperature dependencies of most observables in the 
vicinity of the impurity: details appear in Sections II B 
andlilCl 



2. Finite density of impurities 

The problem of a finite density of impurities is one 
of considerable complexity, and we will only address a 
particular aspect for which we have new, significant, and 
experimentally testable predictions. We will not settle 
the issue of whether Fig ||a or b, or some other possibility, 
is the correct phase diagram. 

We will be interested in dynamical properties of the 
region, s > s c , as characterized by the susceptibility at 
the antiferromagnetic wavevector 



XQ^n) 



d d xd d x' 



dr((j) a (x,T)<j) a (x',0)}e 1 ' 



(1.39) 

where V is the volume of the system, and aj n is a Mat- 
subara imaginary frequency. In the absence of impurities 
in an insulating antifcrromagufit, this susceptibility has 
a pole at the spin gap energyEJ (at T = 0): 



XqH 



A 



A 2 - (huJ + iO+) 2, 



(1.40) 



where u> is now a real frequency, + represents a positive 
infinitesimal, and A is a residue determining the spec- 
t ral w eight of the spin-1 collective mode. The form in 
(1.40) is valid only in the vicinity of the pole, and quite 
complicated structures appear well away from the pole 
frequency. Such a pole will also be present in a d-wave 
superconductor if momentum conservation prohibits de- 
cay of <ft a into low-energy fermionic ^ quasiparitcles. The 
case for which such a decay is allowed is considered in Ap- 
pendix [A| we show there that ( 1.40D is replaced by the 
more general scaling form 



A 



(1.41) 



scale A. Such a pole (or broadened pole) has immedi- 
ate experimental consequences: it leads to a 'resonance 
peak' in the neutron, scattering cross-section, as is seen 
in YBa 2 Cu30 7 EZiyyMEI At current experimental res- 
olution, no intrinsic broadening has been observed at low 
temperatures, and so it is reasonable in a first theory to 
work with the sharp pole in (1.40), as we do in the body 
of this paper. 

Our p rimary interest will be in the fate of the pole (or 
of (1.41)) in the presence of a finite density of impurities, 
n-mip- This is a property at the energy scale of order A, 
and useful results can be obtained without resolving the 
very low energy properties which distinguish the phase 
diagrams of Fig ^. At small n- lmp , the typical interaction 
D in (1.25) is exponentially small in (A/r) 1 / 2 , and we 
will neglect such effects. Instead, at the energy scales we 
are interested in, we use mean-field approach to account 



for a finite 



but include the full dynamics of the 



interaction between the 4> a and a single impurity: the 
finite rii m p will only modify the density of states of the cj) a 
excitations, and this will be determined self-consistently. 

Let us fir st re view the exact scaling arguments. The 
expression ( 1.39 ) defines a bulk observable which cannot 
acquire any of the anomalous boundary dimensions. Con- 
sequently its deformation due to the impurities is fully 
determined only by the energy scale T define d in (1.24), 



and we have the scaling prediction that (1.40) is modified 
to 



XqM 



A_ 

A 2 



huj r 
a"' A 



(1.42) 



where $ is a fully universal function of its arguments 
(we have restricted attention to T = here — for T > 
there is universal dependence on an additional argument, 
T/A). For the case where there is no bulk coupling to 
fermionic quasiparticles at the lowest e nerg ies, we have 
$(w,Q) = 1/(1 - (u + iO + ) 2 ), and so ( |l.42| ) reduces to 
( 1.40 ) for zero impurity density, T — 0; with a bulk cou- 
pling of 4> a to fer mionic quasiparticles, as in Appen dix 
( |1.42| ) reduces to ( |l.4l| ) at T = 0. The scaling form ( |l.42|) 
applies to the dynamic susceptibility for both possibili- 
ties of the phase diagram in Fig g; if Fig ||b is the correct 
phase diagram, then the scaling function $ will have a 
singularity at the critical value of T/A. 

Notice that the impurities induce broadening at an en- 
ergy scale r ~ 1/A. On the other hand, intrinsic broad- 
ening is suppressed by momentum conservation; in the 
special case that Q = (2K, 2K), intrinsic broadening ex- 
ists at the lowest energies, and then it is at most of order 
A. So for small A, the impurity broadening dominates 
in both cases. In experimental comparisons it may be 
useful to use fits with a linewidth written as a sum of 
intrinsic and extrinsic contributions. Further discussion 
of these issues is in Appendix [A|. 

We will obtain results for the form of $ for non-zero 



where $ is a universal scaling function. The form (1.41) T/A in Section IV in a self-consistent non-crossing ap- 



predicts that in such d-wave superconductors the pole 



in (1.40) will be universally broadened on an energy 



proximation, for the case where the coupling to fermionic 
quasiparticles is irrelevant. As we mentioned earlier, one 
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of our important results will be that the T > lineshape 
has an asymmetric shape, with a tail at high frequen- 
cies. This is a significant prediction of our theory, which 
should be testable in future experiments. The experi- 
ment of Ref. ff| has S r = 1/2, T = 5 meV, T/A = 0.125 
and we display our predicted lineshape for this value of 
T/A in Fi g H (results for other values of T/A appear 
in Section |TV| ). The half- width of the line is approxi- 
mately T, and this is in excellent accord with the mea- 
sured lincwidth of 4.25 meV. The sharp lower threshold 
in the spectral weight is an artifact of the non-crossing 
approximation: consideration of Griffths effects (which, 
i n pri nciple, are contained in the exact scaling function in 
(1.42)), or the contribution of the "irrelevant" fermionic 



excitations (which are not part of the scaling limit ( |l.42[) 
will always lead to a small rounding of the lower edge. 
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FIG. 4: The universal lineshape Im<£> as a function of Tuu/ A for 
d — 2, S r = 1/2, and F/A = 0.125, the values corresponding 
to the experiment of Ref. [L5|. This result has been obtained 
in a self-consistent 'non-crossing approximation' analysis of 
5b + Si m p described in Section |lVj 



D. Outline 

We now sketch the contents of the remainder of the 
paper. Readers not interested in details of the deriva- 
tion of the results may wish to skip ahead to Section [v|. 
Section [n| will obtain renormalization group results in an 
expansion in e = 3 — d for the single impurity problem. 
Similar physical results will be obtained in Section III, 
but in a different large-A" or 'non-crossing' approxima- 
tion (NCA). The advantage of the latter approach is that 
it can be readily extended to obtain a self-consistent the- 
ory of the magnon lineshape broadening; this shall be 
presented in Section IV. We review our main results and 
discuss some further experimental issues in Section [v|. 
The appendices contain extensions and further calcula- 
tional details; in particular, Appendix |A| considers the 
case where momentum conservation allows a bulk cou- 
pling to the fermionic quasiparticles of a d-wave super- 



conductor at the lowest energies, while Appendix |b] gen- 
eralizes our results to antifcrromagnets and d-wave su- 
perconductors with incommensurate spin correlations i.e. 

Q^Tr). 



II. EXPANSION IN e = 3 - d 

The central ingredient behind the universal structure 
of all of our results is the analysis of t he re normalization 
group equations presented in Section [I A . This is per- 
formed order-by-order in an expansion in e = 3 — d. The 
results will be used to obtain expressions for a number of 
physical observables in the subsequent subsections: we 
will consider the paramagnetic region, ,s > s c , in Sec- 
tion II B , an d the magnetically ordered phase, s < s c , in 
Section II C. 



All of the computations in this section will consider 
only the single i mpur ity problem. We will therefore use 
the notation in ( 1.30| ) and also denote 



70,r=0 = 7o- 



(2.1) 



Also in the remainder of this paper, we will use units of 
time, temperature, and length in which 

h = k B =c=l. (2.2) 



A. Renormalization group equations 

We will follow the orthodox field theoretic approachEl 
in obtaining the renormalization group equations of Sb + 
Simp'- generate perturbative expansions for various ob- 
servables correlators, multiply them with renormalization 
factors to cancel poles in e, and finally use the indepen- 
dence of the bare theory on the renormalization scale to 
obtain the scaling equations. This approach is rather 
abstract, but has the advantage of allowing explicit cal- 
culations to two-loop order and establishing important 
scaling relations to all orders. A more physically trans- 
parent 'momentum shell' formalism can be used to obtain 
equivalent results at one loop, but several key properties 
are not delineated at this order. 

It is clearly advantageous to have a diagrammatic 
method which allows expansion of the correlators in pow- 
ers of the couplings go and 70. While the standard time- 
ordered perturbation theory can be used for g , the non- 
linear Berry phase term in 5i mp (and the associated unit 
length constraint on n a ) makes the expansion in 70 more 
intricate. Building on earlier work in the context of the 
Kondo problemLJ, we have developed a diagrammatic 
representation of the perturbative expansion in 70 — this 
is described in Appendix 0. The representation works to 
all orders in 70, but does not have the benefit of a Dyson 
theorem or the cancellation of disconnected diagrams; 
nevertheless, all terms at a given order in perturbation 
theory can be rapidly written down. For instance, the 
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diagrams shown in Fig |5| lead to the following lowest or- 
der representation of the two-point correlator of n a in 
the paramagnetic phase: 






(b) 



FIG. 5: Representation of the correlation function in ( |2,; 
the diagrammatic approach of Appendix |c| The open circles 
represent the external sources, the dashed lines are (j> a prop- 
agators, and the filled circles contribute factors of 70. The 
impurity spin is represented by the single full line along which 
i mag inary time runs periodically from to (3. The correlator 
(fO) is the ratio of the diagrams in (a) to those in (b), with 
the latter representing the normalizing partition function Z. 



S 2 (n a (T)n a (0)) =S(S + 1) 



"7o 2 



J dn J dT 2 D{ Tl - r 2 ) + • • • [2 



.3) 



where r > and D(t) is the two-point <j) a propagator at 
x = 



D(t) = (0 a (O,T)^ Q (O,O)) o 



d d k 



(2Tr) d uj 2 +k 2 + m 2 



(2.4) 



Here the subscript in the correlator indicates that it 
is evaluated to zeroth order in go. However, we have in- 
cluded Hartree-Fock renormalizations in the determina- 
tion of the 'mass' m; these, have been computed in earlier 
work in the e expansionEj, and we quote some limiting 
cases: 



A ; s > s c , T< A 

(10e/33) 1/2 7rT ; T ~» \s - s c \" 



(2.5) 



Recall that A was defined earlier as the exact T = spin 
gap of the host antiferromagnet. The crossover function 
between the two limiting results in (|2.5|) is also known for 



small e, but we refer the reader to the original paperEil 
for explicit details. 

The T = limit of ( j2.3|) must be taken with some care: 
the function D(t) in (|2.4| ) is periodic as a function of r 
with period 0, and so any significant contributions for 
small I r I have periodic images in the region with small 
|/3 — t\. Indeed, at T = 0, and at the critical coupling 
s = s c (where m = A = 0), (2.3) reduces to 



S 2 (n a (T)n a (0))=S(S+l) 



9 

"7o 



where 



G?7"l 



Do(t) = 



dn 



1 



dr 2 )D (ti - t 2 ) 



2.6) 



d d k dui e~ iu}T 
{2n) d 2^k 2 +LU 2 



with 



Sd+i 



r(d/2-l) 
4^72 



(2.7) 



(2.8) 



Simple power counting of the integrals in (2.6) shows that 
they lead to poles in e-this is as expe cted from the tree- 
level scaling dimension of 70 in (fLTfj). In the field theo- 
retic RG these poles have to be cancelled by appropriate 
renormalization factors which we will describe shortly. 
Evaluating the integrals in ( |2.6D , and also those in the 
two-loop corrections to ( |2.6| ) described in Appendix |(], 
we obtain 



S 2 K(tK(0))=5(5+1) 



-(loSd+lT-) { — + - 



1 



c(l-e) 



(2.9) 



We have retained only poles in e in the co-efficient of the 
7q term, as that is all that shall be necessary for our 
subsequent analysis. 

We now describe the structure of the renormalization 
constants that are needed to cancel the poles in (2.G) 
and in other observable correlation functions. First, the 
renormalization factors of the host antiferromagnet can 
only depend upon the bulk theory S\, as a single impurity 
cannot make a thermodynamically significant contribu- 
tions. These are well known and described in numerous 
text books and review articles; we will use here the con- 
ventions of Brezin et alEQ They define the renormalized 
field 4>R a and the dimensionlcss coupling constant, g by 



9o 



Z 2 S, 



d+1 



(2.10) 



where /1 is a renormalization momentum scale, Z is the 
wave- function renormalization factor of the field d> nn Za 
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is a coupling constant renormalization, and 

2 



Sd 



r{d/2)(4w) d / 2 ' 



(2.11) 



The explicit expressions for Z and obtained by Brezin 
et aini, to the order we shall need them are 



Z = 1 



Za = 1 



144e 
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Now let us turn to the boundary renormalization fac- 
tors associated with the presence of the impurity spin. 
As in the bulk, we have wavefunction (Z 1 ) and coupling 
constant (Z 7 ) renormalization factors: 



z n/2 



n-Ra 

^ 2 Z„ 



7o 



(zz>s d+1 y/* 



(2.13) 



where 7 is the re norm alized, dimensionless boundary cou- 
pling. In both (2.10) and ( [2.13 ) we have inserted judi- 
cious factors of the wavefunction renormalizations in the 
redefinitions of the coupling constants go and 70-these 
are determined simply by the powers of the fields mul- 
tiplying the couplings in the action. The action has a 
term 7o0 a fJ. a , and so the renormalization of 70 picks up 
one power each of the wavefunction renormalizations of 
(j> a and n Q ; similar considerations hold for go . 

The expression ( |2.15 ) contains two boundary renor- 
malization factors, b ut so far we have evaluated only one 
correlation function, (2.S). To obtain an independent de- 



termination of Z~, we consider the correlation function 



S((j) a (x,T)n a {0)} 



(2.14) 



which will involve the coupling 70 at leading order, as 
shown in Fig ^a. This leading term evaluates at T = 
and s = s r to 



70^(5+1) 



d d k 



th\r 



(27r) d k 2 



(2.15) 



Let us first consider corrections to ( 2.15 ) from interac- 
tions involving only the boundary coupling 70. It is easy 
to see that the diagrams for these have interactions along 
the impurity spin loop that are in one-to-one correspon- 
dence with those appearing in the evaluation of the two- 
point n a correlator; these are shown in Figs ||, |l[ [l2| 
and lead to fl2.9p . However, such diagrams are obviously 
associated with the wavefunction renormalization, Z', of 
n a . So there are no contributions to the coupling con- 
stant renormalization, Z 7 , from the boundary interac- 
tions alone; this leads to the important result 



1 at 



0. 



(2.16) 



Next, we consider interference between the bulk interac- 
tion, go, and the boundary interaction, 70. Now there 




(a) 




(b) 



FIG. 6: Repr esent ation of contributions to the correla- 
tion function in (2.14) in the diagrammatic approach of Ap- 
pendix^ the diagram in (a) is the zeroth order result, while 
(b) contains diagrams of order "fog. The wavy line is the in- 
teraction go and other conventions are as described in Fig [|. 



are contributions which are specifically associated with 
the renormalization of 70 and these are shown in Fig ||b. 
Evaluating the graphs of Fig [|b at T — and s — s c as 
described in Appendix |^, w e obt ain the following contri- 
bution to the correlator in ( |2.14[) 



9olo 



S(S + 1] 



d d h d d k 2 



d d k 



(2ir) d k 2 



S(S- 



(27r) d (2ir) d k\k 2 (k + k x + fc 2 ) 2 



(2.17) 



After determining, by standard methods, the residue of 



the pole in e associated with the integral in (2.17), we 
can write down the coupling constant renormal: 



ization 



1 



j [g(5+l)-l/3] g72 
6e 



(2.18) 



The above result for Z 7 -pL is smaller by a factor of 3 from 
that quoted earlier by usc.1, and this corrects a numerical 
error in our computation. This correction also leads to 
corres ponding minor changes in ( [2.23] ), (2.24), ( |2.26 ) and 
( p~83| ) below. 

We now have all the ingredients to determine the fi nal 
renormalization factor Z' . Using the definitions (2.13) in 
fl2.9| ), along with the result (2.18), and demanding that 
all poles in e cancel, we obtain 



Z' = l 



2 7 2 



2_ 



(2.19) 



Now the RG beta-functions can be obtained by stan- 
dard manipulationsEa For the bulk coupling g we ob- 
tain from the definition (2.10) and the values ( p,12| ), the 
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known result 



</(j 



-eg l 



-eg- 



<i In ^dlnZ 
ding ding 

llg 2 23 5 3 



G 



12 



(2.20) 



Similarly, we can obtain the beta-function for the bound- 
ary coupling 7 defined as 



c?7 
(/// 



(2.21) 



9o ,7o 



Taking the derivative of the second equation in (2.12) 
at fixed bare couplings go, 70, we get 



1 = BM^l 

2 PK1) 9 7 



^)' ln ^ 



18(7) d\nZ' 



dg 2 9 7 

/3(g) & In Z' /3( 5 )dlnZ , /3( 7 ) 



dp 



7 



(2.22) 



Solving Q2.22D using ( |2.12| , pl^j2~19| , p.20D we obtain to the 
order we are working 



/3(7) 



£7 
2 



r 144 



S-7 3 7T 2 



[5(5+1) -1/3]; 

(2.23) 

the first two terms in this betatfunction agree path the 
earlier results of Sm ith and SiEI I and Senguptacj. The 
beta-functions ( 2.2C ) and (2.23) have infra-red stable 
fixed points at 



7 



6e 
IT 

e 

2 4 
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[5(5 + 1) - 1/3] (2.24) 



(1.33). For s > s c in the T — quantum paramagnet, 77' 
determines the impurity static moment, as was asserted 
in ( 1.36| ) and (1.37). These facts are demonstrated ex- 
plicitly in perturbation theory in Appendix O . 

We close this section by mentioning an interesting bi- 
product of our analysis: we will obtain an exact expo- 
nent for some simpler models considered earlier in the 
literature. Consider the fixed point of the beta-functions 
above with g = and 7 = 7* 7^ 0, so that the bulk 
fluctuations of the <f> a are described by a Gaussian the- 
ory. Such a fixed point is unstable in the infra-red to 
the fixed point we have already considered, and so does 
not generically describe impurities in any antiferromag- 
net. Nevertheless, it is an instructive model to study, 
and closely related models appear in mean-field theories 
of quantum spin glassesc3. The fixed poin t value of 7* 2 
can of course be easily obtained from ( 2.23| ) to second or- 
der in e-this value will have corrections at all orders in e, 
and an exact determination does not seem to be possible. 
Even so, we can obtain an exact result for 77'. This is a 
consequence of the identity ( [2.16 ). Using this result, and 
the fact that all g der ivativ es ca n be neglected at g = 0, 
the expressions ([2.22 ) and (2.26) simplify to 



/3(7) 



e 7 / 1 dlnZ 1 
Y V 2 dln 7 
dhiZ' 



(2.27) 



7=7* 



It now follows immediately from the fact that 0(jy*) = 0, 
that 

77' = e (2.28) 
exactly at the g — 0, 7 = 7* fixed point. 



These fixed point values control all of the universal phys- 
ical quantities computed in this paper. In particular, we 
can now immediately obtain the exponents 77 and 77'. The 
value of 77 is, of course, known previously, and is given by 



77 = (i- 



dlnZ 



9=9" 



did) 

242' 



dhxZ 



dg 



9=9* 



(2.25) 



For the new boundary anomalous dimension 77' we have 
dlnZ' 



77 = n- 



dfx 



9=9" ,7=7* 

,ainZ' a . .dlnZ' 



da 



ff=9*,7=7* 



' - <~< ^ + |f[S(S + l)-l/3]J.(2.26) 

This exponent 77' controls the decay of the two-point n a 
correlations at the T = critical point at s — s c , as in 



B. Spin correlations in the paramagnet 

The results in this section will be limited to s > s c , 
i.e., the right half of Fig [| 

We will be interested in numerous different linear re- 
sponse functions of 5b + 5i mp . To discuss these in some 
generality, we extend the actions by coupling them to a 
number of different external fields: in 5b we make the 
transformation 

{d T ^ a ) 2 -> (dr - ie a p- ( H uP {x)(l>- ( ) 2 - H sa (x)$ a , (2.29) 

while we add to 5; mp the coupling 



5 / dTH imPja n a (T) 



(2.30) 



The field H u (x) is an external magnetic field which varies 
as a function of x only on a scale much larger than the 
lattice spacing of the underlying antiferromagnet. In con- 
trast, H s (x) is a staggered magnetic field which couples 
to the antiferromagnetic order parameter — so it oscillates 
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rapidly on the lattice scale, but has an envelope which 
is slowly varying; such fields cannot be imposed directly 
in the laboratory, but associated response function can 
be measured in NMR and neutron scattering measure- 
ments. Finally H unp is the magnetic field at the location 
of the impurity. So a space- independent, uniform mag- 
netic field H applied to the antiferromagnet corresponds 
to H u (x) — H, H s — 0, and H imp (x) = H. We have 
also taken all fields to be time independent for simplicity 
- it is not difficult to extend our approach to dynamic 
response functions with time-dependent fields. 

We can now consider a total of six different response 
function to these external fields. We tabulate below re- 
sults for these response functions obtained in bare pertur- 
bation theory, to lowest non-trivial order. The computa- 
tion is performed by the methods of Appendix ^ and is 
quite straightforward-we do not explicitly show the Feyn- 
man diagrams associated with these results. The results 
will subse quent ly be interpreted using the RG formulated 
in Section II A. 



Xu,u(^, X ) 



Xn,imp{x) 



T 



S 2 InZ 



3 SH ua (x)5H ua (x') 

5 d (x-x>) 2 -^±l±rtTGHx) 

T 5 2 \nZ 
3 5H ua (x)8Hi mp , a 

5(5 + 1) 2 f d d k 1 d d k 2 e i(fci-**)* 
"7o 



3T '"J (2Tr)d(2n)d e kl e k2 
{e k2 coth(e fcl /2T) - e kl coth(e fe2 /2T)) 



•4J 



2T 



Ximp, imp 



Xs,imp( 2 ') 



T 



£ki£k 2 
5 2 InZ 



3 5H ua (x)5H sa (x') 

10Xu,imp(x)G(x') 

T S 2 lnZ 

3 ^-ffimp,a^-^imp,a 
S(S + 1) , o x 

— ^ — (1 - 

T 



S 2 InZ 



3 5H sa (x)5H imp , a 

-2^f^7o 3 G(x)^ 

T 5 2 lnZ 

3 SH sa (x)5H sa (x') 



5(5 + i; 
3T 



1 2 G(x)G(x') 



(2.31) 



where 



G(x) 



(k 2 



2,1/2 



d d k 



(2n) d e\ ' 



T 



d d k /coth(£ fc /2T) 



{2tt)° 



24 



T 



(2.32) 



The prefactors of 1/3 in Q2.31 ) merely perform the rota- 
tional averaging over the directions in spin space. It is 
important to note that we have only included impurity 
related contributions in the above results: the portions of 
the susceptibilities which are a property of the bulk the- 
ory, 5b, alone have been dropped as they were considered 
in earlier work. 

We will now combine these results into various observ- 
able correlators and discuss their scaling properties. 



1. Impurity susceptibility 

First, let us co nsider the impurity susceptibility, Ximp, 
define d in Jl.ll| ). This is related to the expressions in 
(gH) by 



Ximp — Ximp.imp H~ 2 J d ^Xu,imp(*^) 

+ / d d xd d x'x^ n {x,x') 



5(5+1) 
3T 



1 



7^ 



d d k 



1 



2TJ (2ir) d £ 2 smh 2 {e k /2T) 



(2.33) 



An important property of the above expression is that 
as T — > for s > s c , the term of order Jq is exponen- 
tially suppressed and the response reduces to that of a 
free spin 5: this indicates, as expected, that the total 
magnetic moment associated with the impurity is pre- 
cisely 5 in paramagnetic phase. This last fact can also 
be established directly by computing the magnetization 
in the presence of a finite field H at T = 0; this can be 
done to all ord ers in H by the meth ods to be developed 
in Section II C, and the result ( [2.78 ) holds — we will not 
present the details here. 

Next, we examine the behavior of Ximp as a appr oaches 
s c from above. The momentum integral in ( 2.33| ) is ex- 
ponentially convergent, and so there are no poles in e. 
This implies, as expected from the conservation of total 
spin, that Ximp acquires no anomalous dimensions, and 
.33|) can be written in the scaling form 



Ximp 



1 /A 



(2.34) 



with $i mp a universal scaling function. The latter can be 



evaluated from ( 2.33j ) by setting y = 7* , and by using the 
crossover function for m in Ref. |S0j. Let us con side r such 
an evaluation explicitly at s = s c . Then, from ( |2.5| ), m is 
proportional to T, but with the proportionality constant 
of order e 1 ' 2 . It is therefore tempting, to leading orde r 
in e, to simply set m = in the integrand in (2.33). 
However, this leads to trouble — the integrand in (2.33) 
is infrared singular, and behaves like l/(k 2 + m 2 ) 2, at 
low momentum. So we have to keep a finite value of 
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to, and the correction of order in To in (2.33), which is 
superficially of order e, turns out to be of order e 1 / 2 . By 
such a calc ulation, we find that at s — s c , Ximp takes the 



form ( 1.31 ), with 

5(5 + 1) 
w — z 



33e 
40 



1/2 



0(e) 



(2.35) 



Actually, it is possible to also obtain the O(e) contribu- 
tion above without too much additional difficulty. To do 
this, we have to identify only the contributions to Xi mp at 
order Tg and j^go, which are infrared singular enough to 
reduce the expression from superficial order e 2 to order 
e. It turns out that there is only a single graph for Xu,u 
which can accomplish this, and it contributes 



5(5 +1)5 M0 7 2 



d d k 1 



(2.36) 



to Ximp- Evaluatin g (|2.36 ), and also i denti fying the 0(e) 
contribution from Q2.33 ), we correct ( p. 35 ) to 



Ci = 



5(5 + 1) 



33e 

lb" 



1/2 



7e 
T 



o 



(2.37) 



2. Local susceptibility 



Second, consider the response, at or close to the im- 
purity site, to a local magnetic field applied near the 
impurity site. This is usually denoted by xioc and could 
be measured in a muon spin resonance experiment; we 
define it by 



Xloc 



= z'~ 



Xi 



imp, imp ■ 



(2.38) 



imp, imp 



IS 



We have inserted a prefactor of Z'^ 1 because \ 
the two-point correlator of the n a , and this acquires a 
field renormalization factor Z' -1 / 2 in (2.12). It is now 



easy to verify from ( 2.31 



and (2.19) that the poles in e 



do cancel, and that (2.3S) can be written in the form 



Xlo 



-^U, (- 



(2.39) 



Again $i oc is a universal scaling function; it reaches a 
constant value at zero argument and so xioc diver ges as 
T~ 1+r; at ,s = s c . For s > s c , we see from (2.31) that 
Xioc ~ I/? 1 as T — > 0. We can consider this as arising 
from the overlap of the impurity moment with the to- 
tal, fre ely fl uctuating, moment of 5 that was considered 
below (2.33), and so write 



5(5 + 1) 



s > s c . 



cfrp imp 

From the expressions in ( |2.31 ) , we deduce 

To f ddk 1 
2 J {^) d el 



"imp 



= z'- 1 ' 2 



(2.40) 



(2.41) 



It can now be checked that this value for to;, 



"imp agree s 

precisely with that computed from the definition ( 1.36| ) 
- this is shown in Appendix [5], where TO 2 mp evaluates to 
f5s| ). Also this last result, or ( ^39| ), show that ( |l.37|) 
is obeyed with rf = e Alternatively, we can use the 



methods of Section II C to compute the T = impurity 
magnetization in the presence of a finit e field H, to all 
orders in H, and the result agrees with ( |2.41 ). 



3. Knight shift 

Finally, we measure the space-dependent response to 
a uniform magnetic field, H. This can be measured in 
a NMR experiment as a Knight shift, as was done in 
Refs. HJlJ, with the results indicated in Fig |l|. We are 
considering linear response in iJ, and so such results ra#« 
implicitly valid for H <C T. However, experiments&E-a 
are often in a regime where H is of order T, and so it is 
useful to go beyond linear response, and obtain the full 
H/T dependence of the Knight shift. We will show that 
this can be done using some simple arguments in the spin- 
gap regime (s > s c ) with T, H <§; A, but H/T arbitrary: 
from the knowledge of the linear response in H, the entire 
H/T dependence can be reconstructed. In the quantum 
critical region (T ^S> A, s > s c ) we will be satisfied by 
exploring H <^ T in linear response; in the opposite limit, 
the host antifawmnagnet undergoes a phase transition to 
a canted statelSEil induced by the applied field at H ~ T, 
and we wish to avoid such complications here. 

We can identify three important components of a 
Knight shift: (i) K- lmp: the Knight shift of a nucleus at, or 
very close to, the impurity site; (ii) K s (x), the envelope 
of a Knight shift which oscillates rapidly with the orien- 
tation of the antiferromagnetic order parameter, (i.e., as 
cos(Q.x)); and (Hi) K u (x), the uniform component of the 
Knight shift away from the impurity site. We will con- 
sider these three in turn. We measure the Knight shift 
as simply the mean electronic moment at a particular lo- 
cation, and will drop the factor of the electron-nucleus 
hyperfine coupling. 



(a) K h 



We define K- lrav by 



z ,-l/2 



Xi 



+ d d XXu,imp(x) 



H 



z ,-l/2 



5(5 + 1) 
3T 



7o 2 



d d k /coth(e fc /2T) 



(2Tr) d \ 



i4 



4T£ 2 sinh 2 (e fe /2T) 



H. 



(2.42) 



This has a renormalization factor of only Z'" 1 / 2 because 
it is the correlator of n a with the total magnetization, 
and the latter conserved quantity requires no renormal- 
ization. There is an overall factor of H because the mag- 
netization is induced by the external field, and we are 
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considerin g line ar res pons e. As for (2.39), it can be veri- 
fied from ( 2.31 ) and ( 2.1 9| ) that the poles in e cancel, and 
the result is of the form 



-r,'/2 



rpl- V '/2 



(2.43) 



The scaling function $if im reaches a constant value at 
zero argument, and so K- lnip diverges as T~ 1+v I 2 at s — 
s c . 

For s > s c , Ki mp ~ 1/T as T — > 0, and by the analog 



of the arguments associated with ( 2.40| ) we can now write 



lim K imp = 



S{S+1) 
3T 



''imp 



H 



S > S r . 



(2.44) 



the resulting valu e for m\ mp agrees with ( 2.41 ). We 
can extend ( |2.44 ) to beyond linear response in H for 
T, H A by realizing that the important thermal ex- 
citations in such a regime are simply those of the free 
spin in the external field; so the prefactor of the free spin 
susceptibility in (2.44) is simply the component of the 
free spin wavefunction near the impurity. We can expect 
that the same prefactor applies for arbitrary H/T, and 
so the general response is the same prefactor times the 
free spin magnetization in a field H at temperature T; 
this generalizes (2.44) to 



K imp = m imp SB s (SH/T) ; H,T < A, (2.45) 
where B$ (y) is the familiar Brillouin function for spin 5: 



R ( \ 25 + 1 ,U f 2S + 



1 



1 

25 



coth 



note that B$(y) ~ y for sm all y, and B$ (y - 
Naturally, fl2.45| ) reduces to fl2.44| ) for H -> 
H/T, the impurity Knight shift is 5 times the impurity 
moment m 



25 l 

(2.46) 
oo) = 1. 
For large 



''imp • 



(b) K s 

For the staggered Knight shift, we have the expression 
K s (x) = Z~ 1/2 XsAmp(x) + d d x'xu,s{x',x) H. 



(2.47) 

Notice that now we only have the field scale renormaliza- 
tion of the bulk theory, as we are considering a correlator 
of (j) a (x) with the total conserved spin. To the order we 
are working, we can simply set Z = 1, and then the ex- 
pressions in ( [2.31 ) yield 

5(5+i) r a r d d k 

K s {x) = — — — G(x) 



3T 

coth(£ fc /2r) 



7o 



7o 3 



(2tt) 



ATe 2 k smh\e k /2T) 



H.{2A8) 



Upon using ( p. 13] ) and ( 2.1 9| ) it can be verified that the 
above is free of poles in e, and the result is of the form 



u^ v/2 /A 
Ks(x) = T( l d _ ri)/2 H^ K (-,Tx 



(2.49) 



with $if B a universal function, as is expected from 

general scaling arguments. At the approximation 

we have computed things, rj — 0, and K s (x) ~ 
(7J/T)G(x)[Max(T,A)] £ / 2 . 



As for K lmp , we can generalize (2.49) to obtain the non- 
linear response in H in the spin gap regime (H, T <C A ). 
We first identify a local staggered moment as in (2.44) 



t^o sK ' 3T 5 



S > Sr 



(2.50) 



we have identified the staggered moment as proportional 
to the expectation value of the antiferromagnetic order 
in an applied field in the z direction, and the factor of 
1/5 follows because we have absorbed a factor of 5 in 
definin g the free moment of which (<p z ) is an envelope. 



From ( |2.48 ) we obtain immediately 



,(x))=SG(x) 



7o 



7o 3 



d d k 1 



(2tt 



) d 2e\ 



(2.51) 



this expression for (<p z ( x)) can also be obtained by the 
methods of Section [[I C| by computing the staggered mag- 
netization in the presence of a finite field a t T = 0. 
The f actor in the square brackets, after using (2.13) and 
( |2.19| ), is free of poles in e, and evaluates to order A e / 2 
at the fixed point value for 7. Because of the G(x) fac- 
tor, (<f) z (x)) decays exponentially for large x, and we will 
consider its small x behavior shortly. To ob tain the non- 
linear response we now have the analog of (2.45|): 



K s (x) = (cf> z {x))Bs{SH/T) ; i/,T«A. (2.52) 

So at large H/T, the staggered Knight shift simply mea- 
sures the space-dependent antiferromagnetic moment in- 
duced by the impurity spin. 

We examine the behaviors of K s (x) and (cf) z (x)) for 
small X] th ese ar e specified by the operator product ex- 
pansion in ( 1.38| ): 



lim U) z {x)) 

|x|->0 



???.;, 



| x |(d-l+i,-r ) ')/2 



(2.53) 



where mj mp was obtained in (2.44); a similar result holds 
for the Knight shift in linear response in H but all T: 

Provided the value of rj is such that d — I + 77 — ?/>0 
(which we definitely expect), the staggered Knight shift 
will increase as one approaches the impurity, as seen in 
Refs. j3 Ul0| , |l4| . Indeed, all of the results of this subsection 
are qualitatively consistent with the trends in Fig [l]. 

(c) K u 

The remaining uniform Knight shift is given by 



K u {x) 



Xu,imp (x) + Id x'xu,u{x',x) 



H. 



(2.55) 
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Now no renormalization factors are necessary because 
this is a correlator of the conserved spin of the bulk theory 
with the conserved spin of the total theory, and neither 
of them acquire any anomalous dimensions. Evaluation 
of (2.55Q using (2.31) gives 



. , S(S+1) 2 f d d k x d d k 2 e*fa- k °> 
Ku{x) = — — — 7 



3T iv J (2ir) d (2ir) d e kl e k . 
{e k2 coth(e fel /2T) - e kl coth(e fc2 /2T)) 



At the fixed point of the beta-functions, this is of the 
universal scaling form 



A 



(2.57) 



As in the staggered case, we can go beyond linea r re- 
sponse in H for H, T <C A. Then, following (2.50), we 
define a uniform magnetization (L z )(x)} by 



hm = S{S+1){LM) 

t^o uy ' 3T S 



s > s r 



(2.58) 



and an expression for (L z (x)) fol lows from (2.56). The 
non-linear response, generalizing ( [2.52 ) is 



K u (x) = (L 2 (x))B s (SH/T) ; iJ,T«A. (2.59) 

As before, the small x behavior of the scaling f uncti on 
is controlled by the operator product expansion (1.38): 



and 



lim (L z (x)) ~ |d _J /2 , 

\x\^0 \x\ a *> ' z 



(2.60) 



(2.61) 



C. Spin correlations in the Neel state 

We consider impurity properties when spin rotation 
invariance has been broken in the host antifcrromagnet 
for s < s c . We will restrict our attention to T = 0. 

The field 4> a has an average orientation which we take 
in the z direction. To lowest order in go this expectation 
value is 



( — 

9o 



1/2 



(2.62) 



there is no dependence on x at this order. We will con- 
sider corrections to this result in renormalized perturba- 
tion theory. In preparation, we quote some additional 
properties of the host antiferromagnet. We will need, the 
value of the critical coupling to leading order in ScEJ 



5.9o 
6 



2^ 



d d k 



1 



(27r) d k 2 + lu 2 



(2.63) 



We also define a parameter sq measuring the deviation 
from the critical point 



so- 



(2.64) 



Associated with this is a renormalized s and a bulk renor- 
malization factor Zi 



so 



H. (2.56) and to lowest order in g we hav< 



M 

6e' 



(2.65) 



(2.66) 



The ordering in the host antiferromagnet produces an 
effective field on the impurity spin, and so its fluctua- 
tions are anisotropic. In principle, these can be computed 
order-by-order in 70 by the perturbative methods dis- 
cussed in Appendix ^J. However, spin anisotropy leads to 
a proliferation in the number of diagrams, and we found it 
more convenient to use in alternative approach. Indeed, 
the broken spin rotation invariance implies that tradi- 
tional methods developed for ordered magnets can be 
applied here: we used-t| hg | D yson-Maleev representation 
for the impurity spinESHEJ. A potential problem with 
the Dyson-Maleev representation is that it is designed 
to obtain results order-by-order by an expansion in 1/5, 
while all our results so far have been exact as a function 
of S. However, it turns out that, at T — 0, the Feynman 
graph expansion of the Dyson-Maleev representation of 
our problem remains a perturbation theory in go and 70: 
to each order in these parameters, the results are exact 
as a function of S. This is a consequence of there being 
only a small number of Dsyon-Maleev bosons in every 
intermediate state in the perturbation theory, and these 
can couple only through a limited number of non-linear 
interactions. In contrast, at T > 0, the Dyson-Maleev 
perturbation theory sums over an infinite number of in- 
termediate states with an arbitrary number of bosons, 
and the results are then no longer exact as a function 
of S; in this case the methods of Appendix |c| should be 
used, as they remain exact even for T > 0. Our inter- 
est here is only in T — 0, and so we will use the more 
convenient Dyson-Maleev method-it has the advantage 
of using canonical bosons and so permits use of standard 
time-ordered diagrams, the Dyson theorem, and auto- 
matic cancellation of disconnected diagrams. 

Let us introduce the Dyson-Maleev formulations. We 
label the states of the impurity spin by the occupation 
number of a canonical boson, b. The spin-operator is 
given by 



Sn z 

S(n x - in y ) 
S(n x + iny) 



(2.67) 



Notice that the representation does not appear to re- 
spect the Hermiticity of the spin operators; this is be- 
cause a similarity transformation has been performed on 
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the Hilbert space — we refer the reader to the literature 
for more discussion on this point. It is also convenient 
to define a 'circularly' polarized combination of the bulk 
field <j) a , and to sh ift th e longitudinal component from 
the mean value in (2.62): 



1p = {<f>X + i<j)y)/V2 



—6s 
9o 



1/2 



(2.68) 



Now iS; mp in ( 1.1 3| ) takes the form (we only have a single 
impurity at r = and have dropped the sum over r) 



ch 



-70^(5-6+6) 
-70 ( bty + ip*b ■ 



iP*tfbb 
2S 



,(2.69) 



where it is understood that tp z and ip are evaluated at 
x = in the above. Notice that, at zeroth order, the b 
bosons have finite energy gap of Joy/— 6s/3o — it is this 
gap which stabilizes the perturbation theory and permits 
evaluation of results exact in S at each order in 70 and 
go. Fo r completeness, we also quote the bulk action 5b 
in Jl.4| ) in this representation: 

S h = fdtxj dr \ {{d T 4> z f + c 2 (V^) 2 - 2s<f>£) 

+ \drlP\ 2 + C 2 |V^| 2 + V2 + 



go 

"4! 



^(^+# 2 |^| 2 +4|^r 



(2.70) 



The partition function is now an unrestricted functional 
integral over b(r), 4> z (x,t) and ip{x,r) with weight 
exp(-S b - S imp ). 

Computation of the perturbation theory in 70 and go 
in above representation is a completely straightforward 
application of standard methods. There are a fair num- 
ber of non-linear couplings, and so tabulation of all the 
graphs can be tedious. We will be satisfied here by sim- 
ply quoting the results of perturbation theory and then 
providing a scaling interpretation. We consider a few 
different observables in the following subsections. 



1. Static magnetization 

We consider the static magnetization of the impurity 
site, and also the staggered and uniform moments in the 
host antiferromagnet. 

First, the static magnetization of the impurity, TO; mp , 



defined in (1.36). In bare perturbation theory we obtain 



4mp 



d d k 



2^ 



(fc 2 + lu 2 )(-ilu + 7oV C 6i79o) 2 



(2.71) 



We evaluate the integral, perform the su bstitu ti ons t o 
t he r enorm alized couplings and fields in ( [2.10 ), ( 2.13 ) 



( |2T2 



(2.64) and (12.65), use the rcnormalization constants in 



fl2.18) ), ( |2 19| ) and fl2.66| ), and then expand to the 



appropriate order in e. All poles in e cancel and we obtain 
J In ! 



in. 



1 + 7 



' IE 



2 V *g 



(2.72) 

where 75 is the Euler-Mascheroni constant. Substitut- 
ing the fixed point values 7* 2 = e/2, g* — 6e/ll, and 
exponentiating we see that 



-S x 



(2.73) 



By ( 1.37 ) this defines the exponent rfv/2, and is consis- 
tent with the leading order values v = 1/2, rj' = e. 

Next we turn to the spatial dependence of the static 
moment in the host antiferromagnet. There will be a 
staggered contribution to this which oscillates with the 
local orientation of the antiferromagnetic order, given by 
((j) z (x)). Evaluating the latter in perturbation theory we 
obtain 



5(a)) = 70^ 



d d k 



fJO 

2 



2ir) d k 2 
d d k fduj 



2.s 



-6s 

go 



1/2 



(2ir) d J 2ir (k 2 + uj 2 )(k 2 + u 2 - 2s) 



(2.74) 



This is evaluated by the method described below ( 2.71 ) 
again poles in e cancel and we obtain 



(*)> = 



-3s 
4n 2 q 



1/2 



-In 



(-25) 



(2.75) 



The function F\ has the properties 



Fi(y) = 
Um-Fi(y) ~ 

y->0 

lim Fi(y) ~ 

y-+oo 



47T 

3^ 
1/2/ 



1/2 



K 1/2 (y) 



'/y, 



(2.76) 



where K^^-iiy) is a modified Bessel functio n. A fter sub- 
stituting fixed points values for 7 and g in Q2.75j ) , wf-Sg^i 
that for \x\ — > 00, the moment has the bulk behaviorE§l2l 
~ (— s) 13 with exponent (3 = 1/2 — 3e/22, as expected. In 
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the opposing limit, the result is consistent with t hat ex- 
pected from the operator product expansion ( 1.38| ), given 
in ( 2.53D , to leading order in e. 

A second contribution to the host static magnetization 
is uniform on the lattice scale — this is given by the ex- 
pectation value of the magnetization, (L z (x)), defined in 
(1.20). Perturbation theory yields 



(L z (x)) = 7o 2 S 



d d ki d d k 2 i(kl+k2)x f du_ 
(2-n) d (2it) d J 2tt 

— ILU 



{oj 2 + k 2 )(cj 2 + k 2 )(-iu + 7 OV /-6s/ 5o ) 
t|5 frf d fci d d k 2 e l< - kl+k ^ x 
~ J {2n) d (2ir) d (|fci| + |fc 2 |) 
1 



(M +7ov / - 6s /5o)(|fc2| +7ov / ~ 6s /5o) 



(2.77) 



We will leave this integral in the unevaluated form above: 
s uffice to say that (L z (x)) decays to at large |x|, while 
( 2.60 ) holds for small x. 



are available from the authors. These diagrams were eval- 
uated and expressed in te rms of renormalized couplings 
as described below (2.71) using the Mathematica com- 
puter program. The final expressions were free of poles 
in e, and this was a very strong check on the correctness 
of the computations. At the fixed point value for 7 and 
0, the final result was 
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2V^TT1 



1.00936 



44 



tt5(3V22 + 14tt) 7tt 2 5 2 



330 



165 



(2.79) 



Let us express this result in terms of the spin stiffness of 
the ordered state, p s . We expect XimpX to scale as an 
inverse energy, and tho-oarameter which sets the energy 
scale for general d isE2rc3 



2c 



-l i/(d-i) 



115^ 



(2.80) 



tion 



2. Response to a uniform magnetic field 

In the notation introduced at the beginning of Sec- 
II B, a uniform magnetic field corresponds to 
__{x) = H ai Hi mp ^ a — H a , and H s = 0. We have to 
differentiate between the responses parallel and orthogo- 
nal to the bulk order parameter. 

In the direction parallel to the bulk order (the z di- 
rection), the symmetry of rotations about the z axis is 
preserved, and this means that total spin is a good quan- 
tum numbePn Consequently, the total magnetic moment 
is quantized!^! precisely at 5 



T- 



J\nZ 
5H Z 



S. 



(2.78) 



It can be verified that this holds order-by-order in per- 
turbation theory in 70 and go] the sensitive cancellations 
required to make this happen are consequences of gauge 
invariance. This magnetic moment is pinned to the direc- 
tion of the bulk antiferromagnetic order and is not free 
to rotate-so unlike the situation in the paramagnet, it 
does not contribute a Curie susceptibility. Indeed, the 
direction of the bulk antiferromagnetic order is invari- 
ably pinned by very small anisotropics which are always 
present; in such a situation, the total impurity moment 
is also static, and longitudinal susceptibility is zero. 

Now consider the response to a field in a direction (say 
x) transverse to the antiferromagnetic order, \±- Such 
a field induces numerous additional terms in the action 



<5b+< Sjmp w hich can be deduced from ( |2.29D , fl2.30|) , ( |2.67| ) 
and (2.68). We then expanded the partition function to 
second order in H x in a perturbation theory in 70 and go- 
This required a total of 22 one-loop Feynman diagrams. 
We will refrain from listing them here as the computa- 
tions are completely standard — details of these diagrams 



The numerical factors have been chosen for convenience; 
note that in d = 2 p s is simply proportional to p s . The e 
expansion for p s is available in Ref. jsoU86l 



1 



1 



44 



127e 

484 



(2.81) 



If we eliminate s between (2.7E) and (2.81) we see that 
the p dependence also disappears: this verifies that Xi m p_i_ 
and p s are uni versal ly proportional to each other. We 
may generalize ( 1.32 ) to d < 3 by 



C, 



Ximp_L 



l/(d-l)' 
Ps 



(2.82) 



Our results yield the e expansion for the universal con- 
stant C3: 



155 (IIS, 



+1 



i/(d-i) 



1 



1.19299 



22 V 2e 

7r5(3V22 + 14tt) 7tt 2 5 2 
330 + 165 



(2.83) 



III. 



SELF-CONSISTENT NCA ANALYSIS: 
SINGLE IMPURITY 



In this section we complement the RG analysis of Sec- 
tion [n] by a self-consistent diagrammatic approach. This 
new approach will allow us to obtain more detailed dy- 
namic information for the single impurity problem. It 
can also be easily extended to treat the problem with a 
finite density of impurities, and this will be considered in 
Section 0. 
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One way of motivating the analysis is the large- N ap- 
proximation: the symmetry group of the impurity spin 
is extended from SU(2) to SU(JV) and the N -> oo 
limit is taken by a saddle-point of the functional inte- 
gral. Alternatively, the resulting saddle-point equations 
can also be interpreted as the summation of all 'non- 
crossing' Feynman diagrams — this is the so-called non- 



crossing approximationOEj (NCA). While the NCA and 
large-TV approaches are equivalent for the single impu- 
rity problem, this will no longer be the case in the many- 
impurity analysis of Section ^ — there we will use the 
NCA, but will not have 1/N as a control parameter. 

It is worth pointing out that 1 /N is the only small pa- 
rameter in the present single impurity analysis; we are 
not restricted to small e = 3 — d or to a perturbation the- 
ory in the coupling between bulk and impurity. Indeed, 
we will sum diagrams to all orders in the latter coupling, 
and can work directly d = 2. However, we will formulate 
results in arbitrary d to allow for a comparison with the 
expressions obtained in the e expansion of Section 0. 



Hamiltonian formulation 



It is convenient to present the NCA analysis in a 



Hamiltonian formulation of Sb + Su 



The bulk sys- 



-'lmp ■ 

tern shall be represented by a Heisenberg model of spins 
i on a regular two-dimensional lattice : a concr ete ex- 
ample is an array of coupled laddersc3'ESL3Lilt!3c3. At 
zero temperature the bulk system can be driven from a 
paramagnetic state with energy gap A to a Neel state by 
varying a coupling constant s. 

For an explicit derivation we assume that the paramag- 
netic phase of the bulk is dimerized. Its excitations shall 
be described using the bond-operator formalismEa where 
the spins of each pair forming a dimer are represented 
by bosonic singlet (si) and triplet (ii a , a — x,y, z) bond 
operators: 



QOt 



1 



-J. _ 



(3.1) 



Note that the index i here labels bonds, i.e., pairs of lat- 
tice sites, while the subscripts 1, 2 identify the two spins 
in each pair. The bulk Hamiltonian can now be expressed 
terms of the bond operators. Using standard mean- field- 
type approximations (i.e., condensation of singlet opera- 
tors and decoupling of quartic triplet terms) one arrives 
at 



K = J 



E (McAc* + ^ (MJL ka + h.c.) ) ci.2) 

k . a 



which contains only bilinear terms in t. J denotes a char- 
acteristic bulk energy scale (i.e., the nearest- neighbor 
coupling constant). The dimensionless functions A k and 
-Bk contain the geometry of the system and depend upon 
the coupling constant s and the mean-field parameters. 



7i h can be easily diagonalized by a Bogoliubov transfor- 
mation leading to 



Ti-b = E £ k*L*ka + const 



(3.3) 



k , 



where the new triplet bosons are defined through — 
""k^k + v^_ k , Uk and v k are the Bogoliubov coefficients 
vl = ±1/2 + J^4k/2ek, and ek denotes the dis- 



with it k , t>k 



Bl). Note 



persion of the triplet modes, e k = J (A 
that ^4k and £?k are finite, smooth functions of k. For a 
detailed discussion of appropriate mean-field calculations 
see e.g. Refs.||@ 

The spin-1 excitations t^ a in (3.3) appear to be non- 



interacting bosons. However, this is somewhat deceptive: 
interactions between these bosons where necessary to ob- 
tain the self-consistent dispersion ek, and these interac- 
tions also lead to a crucial temperature dependence in 
£k- The dispersion ek in the paramagnetic state is given 



bye k = 



k at small k (the velocity c = 1 in our 



conventions), whe re m denotes the renormalized mass 



introduced in (2.4), (2.5). At zero temperature we have 
to = A, the spin gap. For T > 0, the int era ction go in 
<Sb led to the temperature dependence in (2.5); similarly 
here we find that the solution of the mean-field equations 
leads to the scaling form 



m = 



T 



(3.4) 



where $ m is a universal scalipg function. Expressions for 
$ m in general d are availableEj, and in d = 2 we have the 
simple explicit resultEZl: 



$ m (A) = 2sinlT 



.A/2 



(3.5) 



Note that in the quantum-critical region, T A , m = 
(21n((\/5 + 1)/2))T, a dependence similar to (2.5) in the 
e expansion. 



In analogy to Section IB we now introduce magnetic 
impurities into the bulk system. Here we restrict our- 
selves to a single impurity spin | coupled to the bulk at 
(bond) site with a coupling constant K: 

Wimp = K i?a(4a + ^Oa) 
i 

K 



ka: 



E^(lL+W- (3-6) 



fk 



Here S a are the components of the impurity spin, and 
N s is the number of dimer sites of the bulk lattice. Note 
that we have omitted terms of the form 6a^3"y Sot t pt'y , Such 



terms correspond to the i rrele vant coupling in (1.21). 
The second identity in eq. (3.6) holds for ek <C J. 

At this point it is useful to establish the correspondence 
between the operators in 7ib+^imp and the fields in S D + 
Simp- The triplet bosons tka near the antiferromagnetic 
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wavevector Q = (tt, tt) represent the fluctuations of the 
antifcrromagnetic order parameter; they are related to 
the (p a fields by 

(*k«+4a)W_ 1 



= ^=j d d xe^-^ a (x,T). (3.7) 



The impurity spin operator S a ( r) directly corresponds 
to Sn a (T) appearing in eq. ( |l.l3| ). 



B. Large-iV limit and non-crossing approximation 

To allow for a controlled approximation we generalize 
the impurity spin symmetry to SU(iV) and consider the 
large- N limit. All results of this and the following subsec- 
tions are limited to s > s c , i.e., to a bulk phase without 
spontaneously broken symmetry. The impurity spin is 
represented by auxiliary fermions /„ (v — 1, ...,N), so 
the N 2 — 1 traceless components of the spin S can be 
written as S vlx = /J/^ — S Vfi /2. A chemical potential Ao 
is introduced to enforce the constraint flf v = N/2; 
for N = 2 this corresponds to S = 1/2, and all results in 
this and the next section will be restricted to this value 
of the spin. The generalization of the bulk system to 
SU(iV) symmetry leads to N 2 — 1 triplet operators t vfl . 
The Hamiltonian for the impurity takes the form 



K 



2JA V 



+ Xo (^2 ftU -N/2 



(3.8) 



We are interested in the dynamics of the impurity spin 
arising from the coupling to the bulk. Feedback effects 
are not present for the case of one impurity; they will 
be discussed in Section IV. The imaginary time Green's 
function for the auxiliary fermions f v is introduced as 



G f , v (r) = -<T/„(r)/t(0)> . 



(3.9) 



Here, T is the time ordering operator in imaginary time. 
If we set K = 0, the unperturbed Green's function is 
simply given by Gfl(iu n ) = l/(iw„ - A ) where u> n 
2 a + 1)tt/{3 denotes a fermionic Matsubara frequency 
and P = l/T. Furthermore we introduce the Green's 
functions for the bulk bosons: 



G t (k,r) 



(T*k(r)4(0)) 



(3.10) 



From the bulk Hamiltonian TCb ( |3.3| ) the Fourier trans- 
form of Gt is found to be Gt(k, iv n ) — 1/ (iv n — et) with 
v n = 2nir/ (3 representing a bosonic Matsubara frequency. 

In order to calculate the self-energy S/^ of the / parti- 
cles arising from the interaction with the bulk-bpsons we 
employ a non-crossing approximation (NCA)E3i2a. This 
amounts to the summation of all self-energy diagrams 



with non-crossing boson lines. The NCA approach can 
be derived from a saddle-point principle in the large-iV 
limit, see e.g. Refs. 77 87. Within NCA we obtain the 



following equations for the Green's functions: 

A* 

where the self-energies are defined by: 
Gj^(iu n ) = iuj n - A - Ef ;U (iu> n ) . 



(3.11) 



(3.12) 



The quantity Gu,1oc{t) represents the bulk fluctuations 
at the impurity site, it is given by the bulk Green's func- 
tion of the operator sum (t + P) together with the mo- 
mentum dependence of the vertex in (3.S): 



G 



tt.loc 



(r) 



(3.13) 



G«(k,r) = - 



2JA V 
£k 



(T (t k + il k )(r) (i k +il k )+(0)) 



Here we have suppressed all SU(iV) indices since the bulk 
is assumed to be isotropic (no broken symmetry). We 
remind the reader that Gtt(k,r) and Gu,\oc{t) are related 
by Fo urier t ran sforms to the propagators of the <j) a field in 
( 1.25 ) and (2.4). The Fourier transform of Gtt.\oc follows 
from (|3.3|): 



G 



tt,\oc 



fan) 



-E 



(3.14) 



Finally, Ao is determined by the equation: 



Nq , 



qo= 2 



(3.15) 



The symmetry of the NCA equations permits solutions 
for the Green's functions which are independent of v\ 
therefore we drop the SU(iV) index from now on. (Note 
that this symmetry would not hold in a magnetically or- 
dered phase of the bulk.) 

We introduce spectral densities corresponding to G/ 
and G ttj ioc, 



Pf(uj) = ImG/(o; 

1 



i0 H 



Pit.iocM = ImG tt ,ioc(^ + i0 + ) (3.16) 

7T 

where ptt, loc(^) = — Ptt,\oc{— follows from the defini- 
tion of Gu,\oc(j). From this symmetry property of Gu it 
can be easily shown that for qo = h the chemical poten- 
tial Ao is zero for all temperatures and values of m (see 
appendix ^) . 
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C. Scaling analysis of the NCA equations 

The system of NCA equations can be analyzed in the 
low-temperature, low-energy regime where all energies 
are well below a high-energy cut-off set by the lattice (e.g. 
the nearest- neighbor exchange J) . In the scaling limit the 
momentum sums can be transformed to integrals with an 
upper cut-off A, so we consider A <§; J, T< J, A-^oo. 
For d < 3 there are no ultraviolet divergences, i.e., the 
results are independent of A for large A. 

We start with the bulk spin fluctuations. The triplet 
modes are gapped with an effective mass m« J, their 
dispersion near Q = (n, 7r) is given by e 2 = c 2 fc 2 + m 2 . 
(Note that the momentum k is now measured relative 
to the antiferromagnetic wavevector Q, and we will ex- 
plicitly display the velocity c here and in the following.) 
The scaling lim it of the local bulk Green's function can 
be found from (3.14). Note that the dominant contribu- 
tions to the momentum integral arise near k = 0, so the 
smooth function can be replaced by its k = value 
A . 



Gtt,\oc{i v n) 



2JA S d ir 



(«2 



2N(l-e)/2 



c d sin(7rd/2) 

2JA S 2 , 2 . 2x 

— My n + m ) 



, (3.17) 



(d = 2) 



with Sd defined in (2.11). We have dropped an addi- 
tive contribution to Gtt,\oc, dependent on A, but non- 
singular as a function of frequency. This is made clearer 
by a Fourier transform to the time domain, where the 
A-dependent contributions are negligible at long times: 
At T — 0, the Fourier transform can be written as 



Gu,\oc{t) 



JA S d n 
c d sin(7rd/2) 



m*- 1 $; ( (rm) (3.18) 



where the scaling function &' tt depends on the product 
t = mi only. It has the following asymptotic behavior: 



T < 1 



e -r/ f (3-e)/2 f>>1 



(3.19) 



From this it is clear that the contribution arising from the 
upper cut-off A falls off as e _cAr , and becomes ~ S(t) for 
A — > oo, which will drop out of the scaling equation (see 
below ) for qo — i. For finite temperature, the scaling 
form ( |3.18| ) has to be replaced by 



G 



tt.Xoc 



(r) 



* JA ? S f!, T d ^ 9 U (tT, %) (3.20) 
c d sin(7rd/2) " V ' TJ y ' 



and so as T — > we must have 

T d ~^ tt (tT, ^) = m d -^' tt (rm) . (3.21) 

If we insert the expression for Gtt,\oc into the NCA equa- 
tions it turns out that the dimensionful factors can be 



combined into a single energy scale K which plays a role 
similar to the renormalization scale [i in the RG analysis: 



K 



K 2 J 2A S d n 



c d sin(7rd/2) ' 



K 2 J 



2A S 2 [d = 2) 



(3.22) 



For convenience we have absorbed the dimensionless 
quantities Ao, S d into K, too. The solution of the NCA 
equations can be written in terms of two-parameter scal- 
ing functions: 



O/M - (I 



e/2 



1 



T l-e/2 K e/2 



/ LU 771 \ 

It' TJ 



(3.23) 



Here, $g and $ p are universal scaling functions which do 
not depend on the microscopic details of the bulk o r the 
cutoff. Note that the scaling form for m in (3.4), (3.5) 
can be used to write the scaling functions in terms of the 
arguments lo/T and A/T. The Fourier transform of the 
scaling function $g is given by the solution of the NCA 
equations ( [3.11 ), ( |3.12| ) in the scaling limit: 



^aiiOjn)- 1 = J2 - &n)*Q(iu' n ) (3.24) 



where Q n = uj n /T, and we have used that Ao = 0. For 
shortness we have dropped the argument m/T in the scal- 
ing functions. (For general values of go 7^ \ one finds 
Ao — £/(mjo) — > for T — * 0, i.e., only £/(mj„) — Ej(iwo) 
obeys a scaling form, and any dependence on the cut-off 
A is absorbed in the chemical poten tial Ao.) 

In general, the solution of ( 3.24 ) must be found nu- 
merically. However, some special cases which permit an 
analytical solution will be discussed in turn. 



1. m = 

We first consider the bulk critical point (s = s c ) at zero 
temperature. The triplet modes are gapless and follow 
a linear dispersion = ck. It can be easily shown that 
t he lo ng-time behavior of the local bath Green's function 
( |3.18 ) is given by 



G 



tt.loi 



(3.25) 



By matching the u> powers in the NCA equations it is 
found that the Green's function Gf has the following 
scaling behavior in the zero-temperature limit: 



Gf(r) 



1 



(AV) £ / 2 



e/2-l 



K e/2 



(3.26) 



More generally, these relations are valid provided that 
T< w C A and 1/A « r « 1/T. 
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It is possible to obtain a complete solutionis of the 
NCA equations in the scaling limit for rn = at low 
temperature which means that we take u>,T — > while 
keeping Co — lo/T finite. We note that in general these so- 
lutions do not correspond to a physical situation since, at 
non-zero temperatures and s — s c , bulk interactions al- 
ways lead to a finite effective mass m ~ T. It is, however, 
instructive to display the complete solutions for m = 0. 
By explicitly performing the Fourier transformations and 
using the scaling form 



Gu,\ oc {t) 



2JA S d+1 (irT) d - 
c d [sin(7rrT)] d_1 



(3.27) 



it can be shownEl that the scaling functions have the 
form 



$ p (u>) = A f e*/ 2 



(27rr / 2 -i r (| + f)r(|-f) 



$ G (f) = -Af ( — 
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(3.28) 



with f = tT, uj — uj/T. The amplitude Af is a 
constant depending on d only. Periodicity requires 
Gf(r + (3) = —Gf(r). Note that the above scaling 
function has a conformal-invariant form, i.e., the finite- 
temperature Green's function follows from the T = 
Green's function by applwiiig the conformal transforma- 
tion z = exp(i27rr//3)ca'Eil. This holds for the value 
qo = h\ there the system obeys particle- hole symmetry 
among the pseudo-fermions / under f< «-> /, which fol- 
lows from the fact that the average number of fcrmions 
per "flavor" is \ . The connection between the NCA-type 
approximation and conformal field theory has been dis- 
cussed in Ref. fnj see also Ref. (76[ 



2. m > 0, T = 

For the case of non-zero m we have not succeeded in 
finding a complete analytical solution of the NCA equa- 
tions. Here we shall briefly discuss the general behavior 
of Gf at zero temperature and m — A > 0. The asymp- 
totic expressions for Gtt,\oc ( |3.19| ) allow to determine the 
asymptotic behavior of G/(t): 



Gf(r) 



'A\ f/2 



const t 



T < 1 
> I 



(3.29) 



where now r = rA. In the long-time limit the Green's 
function G/(r) decays to a finite value which scales as 
(A/K) e / 2 (see also next subsection). 

The spectral density of the bulk fluctuations, p tt (ui), 
has a gap of size A. At small frequencies u> <C A the 
impurity spin behaves like a free spin which implies that 



Pf(u>) has a J-peak at u> = with a weight of order 
(A/ky /2 < 1. At small, but finite w, inherits 
the gap from ptt(w), «-e., P/(w) = for < |w| < A. For 
| a; | > A, pj{u>) is non-zero; the singularity at w = A has 
the form p f - (w 2 - A 2 ) (1_e) / 2 [ P/ - 1/| ln(w 2 - A 2 )| for 
d = 2]. Furthermore, Pf(u>) has singularities at all integer 
multiple frequencies of A. For large frequencies, u» A, 
it decays as | e / 2 1 . A numerical result for the scaling 
function $' (o>/A) associated with the fermion spectral 
density at d — 2 is shown in Fig. |7[ 




FIG. 7: The zero-temperature scaling function $p(a>/A) 
associated with the fermion spectral density pf as a function 
of lu/A, calculated for d = 2. The arrow denotes the 5-peak 
at uo — with an amplitude of 0.34. 



D. Spin correlations in the paramagnet 

In this section we consider the response of the sys- 
tem to a n external magnetic field H in analogy to Sec- 
We restrict our attention to a homogeneous 
we set the staggered component H s = and 



tion 
field 



[IB 



neglect any spatial dependence of the uniform component 
H u . Within our large- N generalization the field H cou- 
ples to the ./V 2 — 1 components of the impurity spin as well 
as to the bulk bosons. The field will be described by the 
hermitean matrix H u ^, the Hamiltonians are extended 
according to 

(3.30) 



which can be deduced requiring that a homogeneous field 
H = H u = H lmp couples to a conserved quantity. It 
is easily seen that at large-A the susceptibility of the 
bulk alone is of order N whereas all impurity correc- 
tions which will be considered below are of order unity. 
Without loss of generality we assume in the following 
H v ^ = H (6 Uj i6n t 2 + <5i,,2<W)A *- e -> tne ne ld couples to 
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a traceless combination of two 'colors' which reduces to 
the x component of the spin for the usual SU(2) case. 
The susceptibility of a free spin in the large- N limit is 
easily found as 



Xfrco 



(-'free 



(3.31) 



with Cfree = 1/8; the exact answer for SU(2) is, of course, 

Cfree = 1/4. 



As in Section II B we consider response functions Xu,u? 
Xu,im P , and Ximp,im P . The Feynman diagrams being rel- 
evant in the large-iV limit are shown in Fig. K Note 





wmp,imp 



imp 



y (D 
A-u,U 



(2) *• -~ 



FIG. 8: Feynman diagrams for the susceptibility to a uniform 
magnetic field. The open circles are external sources, the filled 
circles represent the interaction which contributes a factor 
K, the full lines are propagators of the auxiliary fermions / 
representing the impurity spin, and the dashed lines are t 
propagators. 

that there are two contributions to %u,u = 2xu*u + 
where Xu,u has four interaction vertices. At the criti- 
cal point (s = s c ) the four diagrams have the following 
low-temperature behavior: 



Au,U rji 



( 2 ) ~ I 



^u.u 



T 



(3.32) 



Ximp, imp 



T i- 



Ximp,u ~ yi_ e /2 



In the gapped phase at low temperatures, T <C A <C K, 
the susceptibilities are 



y (i) 

Au.u 



A 



1 

A 



(2) 
Au.U 



1 

T 



(3.33) 



A\ e/2 1 



Ximp, imp ~ I ^ J j, , Ximp.u ~ I ^ J j, 

We now proceed by combining the results into observ- 
able quantities. First, we ha ve the impurity susceptibil- 
ity; Ximp, as defined in ( |l . 1 1|) . It is given by 



2Ximp,u + Xu ,u ■ 



(3.34) 



From ( |3.32 ), ( 3.33| ) and m, A <C K in the scaling limit 
it follows that the low-temperature behavior of Ximp is 



dominated by Xu,u, whereas the other contributions can 
be considered as corrections to scaling. The impurity sus- 
ceptibility takes the scaling form quoted in (2.34), with 
the universal scaling function <I>i m p given by 



Y (i) 

Au.u 



A 

imp ( TjT 
1 



T 



$(i) 



u,u 

A 

T 



T 



, x (2) 

' Au,U 



(2) 



I $ (2) (± 

rji U,U I rp 



(3.35) 



Consistent with our predictions, there are no anomalous 
dimensions involved in Ximp, the reduced coupling con- 
stant K and the bulk energy scale J have dropped out. 
The contributions to the susceptibility can be evaluated 
numerically. At the critical coupling, s = s c , we find 
Ximp to have a Curie form , Ximp = Ci/hsT, as predicted 
in the introduction (1.31), with a universal prefactor C± 
defining an effective spin, as discussed in the caption to 
Fig H In d = 2, the large N value of Ci is 

Cl = Cfree [0.59 ±0.01], , (3.36) 



where Cf roe is defined in ( 3.31 ) 



In the gapped phase at T <C A, it is seen from (3.33) 



that the low-T behavior of Xi 
alone. As shown in Appendix 

/2) 



is determined by Xu?u 
the Curie prefactor of 

(2) 



Xu,u equals that of a free spin, i.e., <i>u : i(A/T — > oo) = 
1/8. This implies that in the gapped phase the impu- 
rity contributes to the total uniform susceptibility in the 
same manner as a free spin, i.e., we have shown the con- 
finement of the impurity spin. 

As a second observable, we again consider the local 
susceptibility which is defined as the response of the sys- 
tem to a magnetic field applied to the impurity spin only. 
It is given by xioc = Ximp, imp, and is proportional to the 
autocorrelation function of the of the impurity spin at 
zero frequency. In the large- N limit it is simply given by 



Xioc = - / dTG f (r)G f (-T). (3.37) 

If we evaluate this integral using the scaling solution for 
Gf we arrive at 



Xlo 



It-' /A 



(3.38) 



with the universal scaling function $i oc (A/T) determined 
by $ g (tT, A/T) from Eq. (|3~2^): 



loc 



df $g | r. 



T J J *-* V T 
The asymptotic behavior of $i oc (A/T) is found to be 



(3.39) 



*loc(A) 



const A <C 1 



A f 



A> 1 



(3.40) 



The local impurity susceptibility diverges at s — s c as 
T -1+e , similarly the impurity spin autocorrelation func- 
tion decays as r _e at the critical coupling. This implies 
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that we have rj' — e in the large- N limit. For T <C A, 
Xioc follows a Curie law with a prefactor of (A/K) e which 
corresponds to a remnant impurity moment of 



<-/2 



mV/2 



(3.41) 



consistent with the prediction (1.37) 



IV. 



SELF-CONSISTENT NCA ANALYSIS: 
FINITE IMPURITY DENSITY 



In this section we will discuss the case of a finite 
density of magnetic impurities using the formalism pre- 
sented in Section III. A finite impurity density can lead 
to va riet y of phen omena which were discussed in Sec- 
tions [B and IC2; we will restrict our attention here to 



the scattering of the bulk spin excitations from the im- 
purities. Interactions between the impurities which even- 
tually may lead to an ordered state at very low temper- 
atures will be neglected. We shall assume T C A and 
work directly in d = 2 and real frequency; all Green's 
functions are retarded quantities, lu implicitly carries a 
positive imaginary part everywhere. All the computa- 
tions in this section will be for the case in which all im- 
purities have spin S r = 1/2. 

The important new effect from a finite n; mp is that bulk 
spin excitations described by the Green's function G t ac- 
quire a self-energy from the scattering at the impurities — 
this changes the local density of states of the bulk exci- 
tations. We will determine the modified density of states 
self-consistently in this subsection. 

The lowest-order scattering process for the t bosons 
is described by the T-matrix T t (w) given by the bubble 
diagram of fermion Green's functions: 

T t (u) = K 2 Tj2Gf^n)G f {iv - iu n ) . (4.1) 



The interaction Hi mp gives rise a momentum-dependent 
self-energy for the t bosons. If we neglect higher-order 
scattering contributions (Born approximation) the self- 
energy is given by 



(4.2) 



where n- lmp is the impurity concentration, and we have 
already performed an averaging over the (randomly dis- 
tributed) impurities in the sample. The Green's function 
G t t ( k, uj) in the presence of scattering can be written as 
[cf. ( gig )] 



G«(k,w) 



AJA V 



lo 2 - [ek + Ef(k,w) 



In the scaling limit we can replace Ak by Aq as above 
leading to 



4JA 



(4.4) 



where the self-energy £ tt (<C £fe) is now momentum- 
independent, it is given by 

E«(w) = 2e k Z t {k,Lo) = 4n imp JA T t (u>) . (4.5) 

The momentum integral can be performed in the scaling 
limit 



G 



ttAoc 



2JA S 2 



ln(m 2 + E tt H-w 2 ) (4.6) 



which replaces the expression (3.17) found at zero impu 



rity concentration. We have used ej. = c k + m . It is 
important to note that the effective mass m acquires a 
renormalization from the scattering processes since the 
Hartree contribution to m contains a boson propagator 
which is .affected by the impurities; m has to be deter- 
mined bycJ 



duj d d k t 
2^(2^l G " (fc ' w) 



G$>(k,u>j)=0 (4.7) 



where G$ denotes the Green's f unct ion ( 3.17 ) in the 
absence of impurities. Equations (4.1), ( |4.5D , (4.6) de- 
termine Gu from Gf, on the other hand Gf is connected 



with G tt through the NCA equations Q3.ll] ), Q3.12Q . So- 
lutions have to be found self-consistently, this procedure 
is similar to a self-consistent Born approximation. 

It is possible to write the solutions in terms of scal- 
ing functions, where the only additional paramet er in - 
troduced by the impurities is the energy scale T (1.24) 
defined in Section I. We restrict the presentation to zero 
temperature, and introduce scaling functions as follows: 



Gtt, loc(u) 

G f {u 



2JA S 2 fu r 
" 1 A' A 



) = (A^r 1/2 *' G (p§ 



A 2 $' 



A' A 



(4.8) 



With the notations m = m/A, uj = y/A , and T = T/A 
we can re- write the equations (p~l|), (|4.5|), (4.6): 



Im$' s (^f) = ^ / ^ilm^(o) 



& tt {Q,T) - ln[m + ^(<D,f)-<D 2 ] , (4.9) 

x Im $' G (Qi,f) [n F {u>x) -n F {u>-u>x)] , 
and the NCA equations ( flj] ), (|1|): 

du>i - 



(4.3) Im($' G (w,r) 



-Im$' G (u) -u>i,T) (4.10) 



x Im <5>' tt (wi, T) [n B {ux) + n F (u> - Wi)] . 

The functions n F and ns are Fermi and Bose func- 
tions in the zero-temperature limit, n F (x) = @( — %), 
ub{x) = —Q(—x), where Q(x) is the step function. $ G 
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and $2 are analytic functions of ui in the upper-half 
complex plane, so their real parts are given by Kramers- 
Kronig relations. Finally, the renormal ized mass m is 
determined by the constraint equation (4.7) writte n in 
terms of The system of equations (l.£), (4-10) can 



be solved self-consistently for a given value of r/A, and 
universally describes the spin dynamics of both the host 
antiferromagnet and the impurities. For Y = the equa- 
tions of course redu ce to t he ze ro-temperature limit of 
the NCA equations (|3.1l| ), ( |3.12| ) quoted in Section [ffl. 



We now turn to the broadening of the boson spin-1 
collective mode peak seen e.g. in a neutron scattering 
experiment. Therefore we have to calculate the dynamic 
spin susceptibility which is identical to the Green's func- 
tion Gtt(k, ui). We consider the dynamic susceptibility 
at the antiferromagnet wavevector 



A 



A 2 -cj 2 + £ tt (w) 



(4.11) 



The scaling form for £ tt in (4.8) implies that xq( w ) m ~ 
deed obeys the scaling form (1.42) with the scaling func- 
tion $ given by 



^(w.r)- 1 = i-(cj + io + ) 2 + $' s (w,r) 



(4.12) 



The results show that Im$(u)) has a hard gap at small 
frequencies, with a gap size of ~ 1—1. 6T. The lineshape is 
strongly asymmetric with a tail at high frequencies. The 
half-width of the resulting line is approximately given by 
T, see Figs. | and |. 
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9: Lineshapes Im$ as a function of hcj/A similar to 
9, but for r/A = 0.05, 0.1, 0.15, and 0.2. 



The asymmetry of the line is a consequence of the in- 
terplay between the gap, A, in the bulk antiferromagnet 
and the inelastic scattering from the impurity spin. The 
bulk magnons self-consistently form an 'impurity band' 
whose lower edge extends to an energy of order T be- 
low A. There is no phase space for scattering below this 
band edge, while frequencies above have scattering states 
readily available — the result is an asymmetric line. 



We now apply the present analysis to the broadening 
of the 'resopance peak' observed recently in Zn-doped 
YBa 2 Cu 3 07B. Using the values he = 0.2 eV, A = 40 
meV, and rii mp = 0.005, we obtain T = 5 meV and 
r/A = 0.125. A linewidth of roughly 5 meV is in ex- 
cellent accord with the measured value of 4.25 meV£3. 
The lineshape observed in Ref. [l5] appears to be asym- 
metric in accordance with our analysis; however, future 
experiments should be able to test this prediction in more 
detail. 



V. CONCLUSIONS 

The object of attention in this paper has been a 
'nearly-critical', two-dimensional, host antiferromagnet 
with a spin-gap A. We described this antiferromagnet by 
a universal, continuum, quantum field theory, Si,: this is 
obtained obtained by deforming a 0(3)-symmetric, 2+1 
dimensional, conformally-invariant fixed poin t by its sin- 
gle relevant operator-the 'mass' term, 2 , in (1.4). Then 
we 'punctured' this continuum field theory by quantum 
impurities at a collection of spatial point {r}. We argued 
that each such impurity was completely characterized by 
a single integer/half-odd-integer valued spin, S r , which 



determined the Berry phase of spin precession, see (1.13). 
In particular, for a non-extensive number of impurities, 
there were no new relevant operators, and so the energy 
scale, A, continued to characterize the dynamics of the 
antiferromagnet in the vicinity of the impurities. A non- 
zero density of impurities, 7ii mp , is a relevant perturba- 
tion, but its impact was universally determined entirely 
by the dimensionless number r/A = n\ mp (hc/ A) 2 , and 
was independent of the magnitude of the bare exchange 
between the impurities and the antiferromagnet. Results 
for many physical properties of this antiferromagnet were 
obtained. 

We also showed how the above results could be ex- 
tended to apply to d-wave superconductors. For the 
case where momentum conservation prohibits a coupling 
of bosonic <f> a field to low energy fermionic Bogoliubov 
quasiparticles, no changes are necessary: the theory for 
the insulating antiferromagnet applies quantitatively to 
the collective spin excitations of the superconductor. The 
situation where there is an efficient coupling of <p a to the 
Bogoliubov quasiparticles is considered in Appendix |A|: 
we show there that the structure of all the scaling forms 
remains unchanged, but there are quantitative changes 
to the numerical values of the scaling functions. 

As a guide to the reader who read Section I and skipped 
the body of the paper, we itemize our main results. For 
the case of a single impurity, our main results were as 
follows. 

• The response to a uniform magnetic field is sum- 
marized i n Fig |3|. The e expansion for the constant 
C\ is in (2.35): it is evident that the series can- 
not be used directly for a quantitative estimate, 
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and we have not attempted any resumm ation anal- 
ysis. The large- N prediction for C\ is in (HH) with 
Cfrcc = 1/4- On the orde red s ide, the e expansion 
for the constant C3 is in (2.83). 



the quantum-critical region is playing a role. In the lat- 



• Local response functions are governed by the known 
bulk anomalous dimension, 77, and the new bound- 
ary anomalous dime nsion rf . The e expansion for 
the latter is in ( 2.26 ). 



The local susceptibility, xioo measures the impu- 
rity site response to a field applied at the impurity 
site only. The e expansion for this i s in (2.38 J2.39 ), 
while the large- N result is in (3.38). 



• The Knight shift measures the localized response to 
a uniform magnetic field applied everywhere. The 
Knight sh ift very close to the impurity site is in 
( 2.42 2.45 ). Away from the impurity site, there 
is a response which oscillates with the orientation 
of the Neel order parameter, as in Fig. 0— the 
envelope o f this oscillating component is given in 
( 2.47 2.49| ). There is also a smaller uniform (ferro- 
ma gnetic) co mponent to the Knight shift — this is 
in ^5^ , ^571) . 



• In the ordered Neel state, t he spontan eous moment 
on the impurity site obeys (1.37,2.72). Away from 
the impu rity site, there is a staggered component 
obeying (2.75), and a uniform component obeying 
07b. 



For the system with a finite density of impurities, we 
obtained results for the magnon lineshape in th e spi n 
gap phase. This line obeys the exact scaling form ( 1.42| 
and NCA resu lts for the sc aling function are in Figs 
and|, and in 

The rest of this section discusses some remaining ex- 
perimental issues. 

One interesting observation of the neutron scattering 
experiments of Ref ^ was that Zn doping caused the 
resonance peak to appear at larger temperatures signif- 
icantly above the superconducting transition tempera- 
ture. We do not have a theory of the superconducting 
order in this paper, and so cannot directly address this 
question. However, this experimental result supports the 
interpretation that the resonance is very much like the 
gapped, spin-1 excitation in the paramagnetic phase of 
an antiferromagnetE3E3. In the "swiss-cheese" madel of 
Zn doping developed by Uemura and collaborators^, such 
antiferromagnetic correlations are enhanced in the vicin- 
ity of the impurity (as in Fig [l] and Refs. ||,[l0 14), and 
should lead to increased stability of the resonance mode 
at higher temperatures. 

Recent measurements of the effective moment on the 
Zn site in YBa 2 Cu306.7 at fairly large temperatures yield 
a value p e g ~ 1.0 which is smaller than the value ex- 
pected for g = 2, S = 1/2, p cS = g- s /S{S + T) = 1.73. 
While some of this reduction is surely due to inter- 
impurity interactions, we speculate that the crossover to 



t er re gime, we expect p e g = g 
( |3.36 ) we obtain p e g « 1.33. 

Next, we discuss the 
work to anal ys e s .. of . . imp -a 
superconductorsE#SE9SBBE 
cent tunneling expcrimcntlij. The former have all been 
carr ied out us ing actions closely related to 5* + S^imp 



3Ci; using the value in 

relationship of our 
.effects in d-wave 
iiHe3 and a re- 



(in (1.26,1.28). However, we have argued in this paper 



that the coupling of the impurity spin to the collective, 
spin-1, mode <p a is a more relevant perturbation. One 
consequence of such a coupling is that the impurity 
spin has structure in its spectral function at frequencies 
whic h arc i nteger multiples of A, as we found in Sec- 
It is therefore natural to expect that some 
structure will appear also in the quasiparticle 
it is interesting that a sideband 



tion 



IIIC2 



of this 

tunneling spectrum: 
at a frequency of order A does appear to be present in 
the tunneling spectrum at the impurity site in Ref. [ll| 
As we move away from the impurity site, the coupling 
between the quasiparticles and the <j) a mode becomes 
negligible, and the frequency A should not be visible 
in tunneling — this also appears to be consistent with 
Ref. |l6|. We note that this sideband phenomenon is 
similar iru-spirit to arguments made by Abanov and 
Chubukovca for photoemission experimental. Also, 
some related computations on the relationship of tunnel- 
ing spectra to antiferromagnetic collective modes appear 
in a recent work by Chubukov et aZo 

In closing, we mention some consequences of our the- 
ory to Ni doping of two-dimensional spin-gap compounds 
of spin-1/2 Cu ions. Each Ni ion has spin-1, and so the 
net uncompensated moment in its vicinity will still have 
S r = 1/2. Our theory, therefore implies that the damping 
of the bulk <f> a quasiparticle mode due to the Ni impuri- 
ties will be the same as that due to non-magnetic Zn im- 
purities, because both depend only the bulk parameters 
and the values of S r ; all differences arise only from ir- 
relevant couplings and these are suppressed by factors of 
A/ J. However, the atomic scale, staggered spin polariza- 
tion will be very different in the two cases; in particular, a 
large contribution to the net moment of S r = 1/2 comes 
from the magnetic moment on the Ni ion itself, while no 
magnetic moment is present on the Zn ions: this implies 
strong differences in the STM spectra. 
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APPENDIX A: FERMIONIC QUASIPARTICLES 
IN d-WAVE SUPERCONDUCTORS 



T his ap pendix will continue the discussion of Sec- 
tion 1 B 1 on the consequences of the gapless, fermionic, 

As 



Bogoliubov excitatio ns in a d-wave superconductor, 
was noted in Section I B 1 , for quasiparticles with vanish- 



ing energy at wavevectors (±K,±K), the key issue was 
whether the antiferromagnetic wavevector, Q, equaled 
(2K,2K) or not. For Q ^ (2K,2K), the bulk cou- 
pling between the antiferromagnetic order parameter, <fi a 
and the fermionic quasiparticles can be neglected at low 
enough energies, and this was assumed in the body of 
the paper. In this appendix, we describe the changes 
necessary for Q = [2K, 2K). 



1. Host superconductor 

The critical theory of an antiferromagnetic order- 
ing transition in a d-wave superconductor with Q = 
(2K, 2Kyj was described recently by Balents, Fisher and 
NayakLa Here we will review their results in a field- 
theoretic approach designed for efficient treatment of the 
impurity problem. 

First, we specify the fermionic action ( 1.26| ) more pre- 
cisely. Let the components of electron annihilation oper- 
ator in the vicinity of the wavevectors (K, K), (-K, K), 
(-K,-K), {K,-K) be f la , ha, ha, ha respectively, 
where a =j, J. is the electron spin component. Now intro- 
duce the 4-component Nambu spinons 'J'i = (fi a , Sab/^) 

and #2 = (ha, Eabflt,) where e a b is an a ntisy mmetric 
tensor with e-fj = 1. Then the full form of (1.26) is 



+ / ^Tj2*U-^ n + k y T* + k x T*)*2- (Al) 



Here r a are Pauli matrices which act in the fermionic 
particle-hole space, the wavevectors k x y have been ro- 
tated by 45 degrees from the axes of the square lat- 
tice, and we have scaled the fermionic velocities to unity 
(following Ref. |5^, we assume that the velocities in the 
two directions are equal near the critical point). With 
Q = (2K,2K), the allowed bulk coupling between the 
antiferromagnetic order parameter <f> a and isEj 



d xdr 



2 ' 



t ^ 1 a y a a T v ^ 1 + H.c. + (1 -> 2] 



(A2) 

where a a are Pauli matrices in spin space. 

The antiferromagnetic ordering transition in the d- 
wave supercondu ctor is described by the action Sb + <S* + 
specified in Jl.4| ), (Al), (A2). We perform its renor- 
malization group analysis in d = 3 — e dimensions by the 



renormalizations (2.10) and 

#1,2 = 

A = 



At Wi2J{ 



/2 



/2' 



(A3) 



At one loop, the results of Ref [5^ translate to the follow- 
ing renor maliz ation constants in the minimal subtraction 
scheme ((2.12) no longer applies): 



z = 


1 - 


4A 2 
e 


5g 2 
144e 


Zi = 


1 + 


Hfl 
6e 


48A 4 


Zf = 


1 - 


3A 2 
2e 




Zx = 


1 - 


A 2 
e 





(A4) 



We have taken the liberty of adding a single additional 
two-loop contribution abov e: th e order g 2 term for Z 
which was also present in (2.12). These translate into 
the beta functions (replacing ( 2.20 )) 



0(a) = -tg + 

/3(A) : 



eA 
Y 



lip 2 
6 

5A^ 
2 



+ 8X 2 g - 48A 4 



(A5) 



The stable fixed point of these equations is 



9 

A* 2 



55 
e 

5' 



(A6) 



and the value of g * rep laces that in ( 2.24 ). The anoma- 
lous dimension in (2.25) is replaced by 



V 



4e 

IT 
3e 

10 



(A7) 



where ry/ is anomalous dimension of the fermions #1,2. 

These fixed point couplings can be used to compute 
the universal correlations of the antiferromagnetic order 
parameter. In particular, in the paramagnetic phase, the 
dynamic s uscep tibility at the antiferromagnetic wavevec- 
tor obeys ( 1.41 ), where A is an energy scale measuring 
the distance from the critical point to long-range antifer- 
romagnetic order. To leading order in e we find 



1 



a 2 -w 2 -s q h' 

with a self-energy, Eq , with imaginary part 

2ttA 2 

Im£Q(w = A) = e — - — . 



(A8) 



(A9) 
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This is the main difference between the antiferromagnetic 
transition in an insulator and a d-wave superconductor 
with Q = (2K, 2K) : the scaling structure in the two 
cases is identical, but the latter has damping in the <p a 
collective mode of the paramagnetic phase even at T = 0. 



2. Quantum impurities 



to 



The 



As discussed in Section IB, 
purities to the bulk theory 

fcrmionic couplings in <S*i mp (Eqn. (IL28I)) continue to 
be irrelevant, and so we discuss the theory with action 



we couple q uant um im 
of Ap pendi x A 1 



Sh+Ssnp+Sv+Svi specified in 0, ( [Til ), Q, flA2|). 



T he a nalysis is essentially identical to that in Sec- 
tion [I A. At the one loop level, there are no new Feyn- 
man diagrams, and only those in Fig ^appear. Moving to 



11 



the two loop level, we find the diagrams in Fig |^, 
and also the single new category of diagrams in Fig 
Evaluating these new diagrams, we find that at this or 




FIG. 10: Additional contribution to the correlator of Fig |^ 
in the theory of Appendix Each propagator of (f> a (repre- 
sented by the dashed line) in Fig ^ acquires a fermion loop 
(represented by the dash-dot line) self energy correction; the 
filled square is the Ao coupling. Only the correction to one 
diagram of Fig ^ is shown here; corresponding corrections to 
the others are implied. 

der, the value of the boundary ren orma lization constant 
Z-y remains unchanged to that in (2.18), while the con- 
stant Z' is modified from (|2.19|) to 



Z' = 1 



e 



7 



2 \2 



+ — + 2 7 Z A 



2 

72 



1 + 2C-In4 



e \ e e 

(A10) 

where C — 0.5772 ... is Euler's cons tant. The renor- 
malization constants in ( 2.18] , A4 A10|) mod ify t he beta 
function for the boundary coupling 7 from ( 2.21 ) to 



0(7) 



.!2 + 7 a + 2A27 _ 7 5 + ^7 

2 ' ' '144 



+^[S(S- 



3 

-i9A 4 7- 



1) - 1/3] - 2AV(1 + 2C- In 4) 

(AH) 



Here i? is a numerical coefficient whose determination 
requires the ord er A term in the bulk renormalization 
factor Z in (A4); the latter has not yet been computed. 
Eqn (All) modifies the fixed point value 7* from (2.24) 



7* 2 = - 
' 10 



(A12) 



Armed with the fixed point values in ( A6| , Al2| ), we can 
now extend all the computations of Section [n] to the the- 
ory of this appendix. There are no changes in the struc- 
ture of the exact scaling forms, but the scaling functions 
and exponents have to be recomputed in the present the- 
ory. In particular, the value of the boundary anomalous 
dimension rf changes from (2.26) to 



(A13) 



APPENDIX B: INCOMMENSURATE SPIN 
CORRELATIONS 

The body of the paper considered antiferromagnets 
with a collinear spin polarization commensurate with the 
underlying lattice. Here we discuss the case of incommen- 
surate spin correlations. 

Several classes of incommensurate magnets need to be 
distinguished. If the incommensurate order is spiral, then 
the spin polarizatiou-is non-collinear, and a very different 
theory is necessaryE3-the quantum paramagnetic phases 
of such magnets possibly have deconfined excitations, 
and their coupling to the impurity spin will be quite dif- 
ferent. We will not consider such a situation here. 

Restricting attention to incommensurate collinear an- 
tiferromagnets, we have to further distinguish two cases. 
Let us describe the spin polarization in the ordered phase 
of such a magnet by 



(S a (x))=Re(N a e i ^ x ), 



(Bl) 



where Q is the incommensurate ordering wavevector, 
and N a is a complex three-component field. Associ- 
ated with any_auch spin polarization, simple symmetry 
considerationsEJ imply that there must be a correspond- 
ing charge density modulation 



(5p(x)) oc cos(2Q -x + tp), 



(B2) 



where (p is an arbitrary phase offset. The final situation 
depends upon the status of the charge density modula- 
tion at the quantum transition at which the magnetic 
order disappears. If the magnetic order (Bl) vanishes 
while the charge modulation in (B2) is preserved on both 
sides of the transition, then the theory in the body of the 
paper continues to apply to the spin excitations. This 
is because we may view the spin modulation as (in a 
sense) commensurate with the underlying charge density 
wave — in other words, the value of tp pins the phase of N a 
and the order parameter becomes effectively real. Only 
in the case where the modulations in (Bl) and (B2) van- 



ish at the same quantum critical point is a new theory 
necessary (the remaining case, in which the order (p32|) 
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vanishes, while (Bl) is preserved is not possible on gen- 
eral symmetry grounds). We have to extend the theory 
5b + 5j mp from a real field n a to a complex field N a . 
We will not enter into this in any detail, but mention the 
changes necessary. In 5b, in addition to the usual quartic 
term, (N*N a ) 2 , a second quartic term, \N a N a \ 2 is per- 
mitted; in (Simp, the coupling 70 becomes complex. The 
phase of "for is a renormalization group invariant (rather 
than just its sign), and its initial value plays a role similar 
to that of a r . 



APPENDIX C: DIAGRAMMATIC 
PERTURBATION THEORY 

We describe a diagrammatic approach to expanding 
arbitrary correlators 5b + 5 mp in powers of 70 and go- 
We will restrict ourselves to a single impurity here, with 
the notations ( 1.30| ) and (2.1). The generalization to the 
many impurity case is straightforward. 

The simplest way to generate this is to use a Hamil- 
tonian representation of the quantum spin in terms of 
Heisenberg spin operators; so we replace 



5? 



S a 7 



(CI) 



where the S a satisfy the standard spin commutation re- 
lations 



[S a , Sp] 



(C2) 



Now we represent the partition function as a time- 
ordered exponential of the interaction terms, expand the 
exponential in powers of 70 and go, and take the trace 
over the Hilbert space of the spin operators. Closely re- 
lated methods .have been developed for the traditional 
Kondo problemEH 

We can reduce the results of such an expansion to a 
simple set of diagrammatic rules. It is important to re- 
member that there is no automatic cancellation of discon- 
nected diagrams in such an approach. So when evaluat- 
ing the expectation value (O) of an arbitrary observable 
O, 



(O) 



Tt(OZ) 
TtZ 



(C3) 



we must separately expand the numerator and denomi- 
nator consistently to each order in 70 and go- The dia- 
grammatic rules can be applied to both these expansions. 

We represent the time evolution of the spin operator 
by a single, directed loop along which imaginary time 
runs periodically from to (3. This loop is shown as a 
full line, and there is one, and only one, loop in each 
diagram. Each factor of 70 is associated with a filled 
circle on the loop at some time tj . Associated with each 
such filled circle we have: 

• A propagator for the (f> a field emerges from the 
filled circle (represented by a dashed line), and 



these are contracted together in powers of go by 
standard 0^ perturbation theory. The go interac- 
tion is represented by a wavy line. 

• A factor of the spin operator S a — these are placed 
in order of the increasing Tj and then the trace is 
taken over the spin Hilbert space. There are also 
factors of S a associated with external sources of 
Sn a in the operator O — these are represented by 
open circles. 

• All the times 73 are integrated over the maximum 
possible interval between to subject to the con- 
straint that they maintain their time ordering along 
the directed loop. 

All of the diagrams contributing at order 7 2 to the nu- 
merator and denominator of the two-point correlator of 
the n a are shown in Fig ||. As an application of the rules 
above, we see that the contribution of the last diagram 
in the numerator (the one with the time labels shown) is 



7o" 



^Tr (SaSpSjp) J dnj dr 2 D( Tl -r 2 ). (C4) 
The trace over spin operators evaluates to 
Tr (sJpSaSp) = (25 + 1)S(S + 1) [S(S + 1) - 1] (C5) 



Other diagrams can be evaluated similarly, and we obtain 
( |2.3D after the quotient has been taken. 

Proceeding to order 70 f° r the same correlator, the 
graphs shown in Fig |ll| are added to the numerator, 
while those shown in Fig [l2| add to the the denominator. 
The new non-trivial spin traces which arise are 



Tr ^SaSpSrySpSaSry^J (25 + l) 

xS(S + l)[S 2 {S + l) 2 - 25(5 + 1) + 1] 

Tr (S a SpSjS a SpSj\ — (25 + 1) 

x5(5 + 1)[5 2 (5 + l) 2 - 35(5 + 1) + 2].(C6) 

We now take the quotient of the additional contribu- 
tions from Fig [IT] and 12|. Although this can be a labo- 
rious process, the work can be shortened by observing 
many cancellations diagrammatically. We have collected 
the diagrams in Fig |ll| in groups of 5 to aid this process; 
all the diagrams within each group differ only by a time 
displacement on the 70 interactions through the external 
vertices. Indeed, the union of the domains of the time in- 
tegrals in each group precisely equals the domain of the 
time integral of a single diagram in Fig [l2[ Were it not 
for the differences between the traces of the spin oper- 
ators in (C6), these contributions to the numerator and 
the denominator would precisely cancel. Consequently, 
we need only keep track of such differences, and can then 
quickly write down the final quotient. The result is the 
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FIG. 11: Contributions to the numerator of the correlator 
(^.3|^at order 70 . These have to be added to the diagrams in 
Fig ga. They have been collected in groups of 5 for reasons 
described in Appendix ^. 



following contribution to the two-point n a correlator at 
order 7q: 



j$S(S + i)J dn J P dT 2 D{ n - r 2 ) 



dr 3 / dr 4 + / dr 3 dr^ + 2 / dr 3 I dr^ 

ii Jti Jt Jt 2 J a Jt 

dr 3 / dr 4 



^(t 3 - n) 



(C7) 



After carefully taking the T = limit of the above ex- 
pression as described above ( |2.6| ), and evaluating the in- 
tegrals using (2.7), we obtain the two-loop contribution 




FIG. 12: Contributions to the denominator of the correlator 
(p.3|)^at order 7$ . These have to be added to the diagrams in 



in (||). 



APPENDIX D: DYNAMIC CORRELATIONS IN 
THE PARAMAGNET 

At T = 0, and s > s c , the two-point correlations of t he 
renormalized n Ra field can be written down from (2.3): 



S 2 (n Ra (T)n Ra (0))=S(S + l) 



-2 7 2 / d Tl / dT 2 D{ Tl - r 2 ) 

JO J -00 



/z\m) 



where have retaine d ter ms accurate only to order 7 2 , and 
at T = 0, D(t) in @) evaluates to 



0(3-d)/2Ad-l q 

J(T)= r(( d -l)/2) (A^V ^-^^), (D2) 

with jffj_ D/ 2 (AT) a m odified Bessel function. The A 
limit of Q is (§r 



We first evaluate (Dl) at the c ritica l point A = 0. 
Performing the integral and using ( |2.19| ) we have 

S 2 (n Ra (r) ■ n Ra (0)) = S(S + 1) (l - 2 7 2 [1 + ln( M r)]) 

(D3) 

to leading order in e- the poles in e have cancelled. For 
large r we can exponentiate this at the fixed point cou- 
pling value 7* 2 = e/2 and find that 



5 2 (n Ka (r)-n fla (0)) =5(5 + 1) 



(D4) 



which predicts 77' = e to leading order in e. 

We now consider the t — > 00 limit of the correlations 
for A > 0: 



lim S 2 (n a (r) ■ n a (0)) = S{S + 1) 

T— >CO 



-2 7o 2 / / dT 2 D( Tl - t 2 ) 

JO J -00 



(D5) 



Introducing new dimensionless coordinates y = A(ri 
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t 2 ), x = A(ri + r 2 ), we write 

lim S 2 {n a {r) ■ n a (0)) = S(S + 1) 

T — *00 

- ^ dyyD(y/A) 



(D6) 



The relevant integral to be evaluated is (with a = (d ■ 
l)/2) 





2 a-l 



(I) 



dyy^K^y) 



1/2 



r(a) 



2 - 2a 



1/2 



r(3/2-a,l) 



(D7) 



In the second line above we have used the asymptotic 
forms of the function K a (y) for very small and large val- 
ues of its arguments, and T(3/2 — a, 1) is the incomplete 
Gamma function. Now putting it all back together 



lim S 2 (n a (T) ■ n a {0)) = S{S + 1) 

T — >00 

-*G(£)' + © , '" r(1A1) ) 

Recalling the definition of Z' ( |2.19| ), we conclude that 

1 



(D8) 



lim S 2 (n Ra (T) ■ n Ra (0)} = S(S + 1) 

T— »00 



(D9) 



Near the fixed point we use 7* 2 = e/2 and exponenti- 
ate to obtai n from (1.36) that mi mp ~ (A//i) e / 2 , which 
agrees with ( |l .37 ) and rj — e. 



APPENDIX E: DETAILS OF THE LARGE- iV 
CALCULATION 

This appendix provides some details of the large-JV 
calculations presented in Section III. 



First we prove that the chemical potential Ao = in 
the particle- hole symmetric case qo = |. From G ti (r) = 
Gu{t)* = G u (-t) and G f (r) = G/(/3 - r) = -G/(-t) 
(for qo — i) follows that Gf{iuj n ) is purely imaginary, 
Gf(ioj n ) = —Gf{iuj n )* = — G f(—iu) n ), and the self- 
energy obeys £/(zu;o) + £/(— mjo) = 0. Therefore the 
chemical potential Ao in the scaling limit is zero for any 
temperature T and effective mass m. 

Second we consider the impurity contributions Xu,u to 
the uniform susceptibility. The diagrams shown in Fig. 



H evaluate to: 
K 2 



x (1) 

Au,U 



y(2) 
Au.U 



2/3 2 



K 

2g5 E Gj(iuj n )G f (iuj n + iv n ) 



E 

4 



Gf(iu> n )Gf(iui n + w n 



x G f (iu> n +iv n )Gf\iv n )Gf\iv n ) . (El) 

The prefactor | arises from the summation over the 
SU(JV) indices together with the form of the field ma- 
trix H Vil defined below ( 3.30[ ). Gj (iv n ) and G\ {iv n ) 
are abbreviations for the product of two and three boson 
Green's functions respectively, together with the momen- 
tum dependence of the vertex; note that the external fre- 
quency is zero here. Due to the antisymmetric structure 
of the field coupling (3.30) both G t {iv n ) and Gt{—iv n ) 
enter the diagrams, with the result being 



GfXiVn) 



G?\iv n ) 



^ 2J lk [G t (k,z^) 2 -G t (k,-^ n ) 2 ] , 



— V 



Ck 



[G t (k,zi/„) 3 + G t (k, 



Evaluation of the momentum integrals using the expres- 
sion for Gt in the scaling limit and performing the Mat- 
subara su mmat io ns in (El) one obtains the behavior 
quoted in flO^) , ( |3~33l) . 

In the gapped phase, T -C A, we expect from the 
confinement of the impurity spin that Ximp = Xfree = 
1/(4/cbT) (with the field conventions used in Sec. Ill Dj ). 
In the low-temperature limit Ximp is dominated by Xu,u, 
since Xa,u remains finite for T <C A, and both Xu,imp an d 
Ximp.imp show a weaker divergence with T — > 0. We prove 

(2) 

Xu,u = Xfree f° r T <C A. We start with the observation 
that 



G?\iv n ) 
Gf\iu n ) 



div u 

-(- 

2 V diu 



G u {iv n 



G t t(if n ) 



(E2) 



which holds at any temperature and independent of the 
dispersion of the bulk bosons. It follows 



K 2 



Gf(iuj n + iv n )G\'(iv n ) 



K 2 



E 



Gf(iuj n + iv n )G t t(w n ) 



(E3) 



From the NCA equations ( 3.11 ), ( 3.12 ) it is easily seen 
that the last term in the [...] brackets equals l/Gf(ioj n ) 

(2) 

in the scaling limit. With this we can express Xa,u hi the 
scaling limit as 



Y (2) = 

Au.U 



^E G K-«) 



diuj. 



G f 1 {iu n ) 



(E4) 
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The dominant contributions to the frequency sum (for 
T — > 0) arise from low frequencies u> n ~ T <C A. 
Therefore it is safe to use G/ ~ l/iu) n (which holds 
for small tu n provided that T <C A). In turn we have 
(d/diw ra )G7 1 (iw„) = (l/iw n )G7 (mj„) which implies 
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